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^S) , The Partition function of two Hard Spheres in a Hard Wall Pore is studied appealing to a graph 

representation. The exact evaluation of the canonical partition function, and the one-body distribu- 
tion function, in three different shaped pores are achieved. The analyzed simple geometries are the 
cuboidal, cylindrical and ellipsoidal cavities. Results have been compared with two previously stud- 
ied geometries, the spherical pore and the spherical pore with a hard core. The search of common 
features in the analytic structure of the partition functions in terms of their length parameters and 
their volumes, surface area, edges length and curvatures is addressed too. A general framework for 
the exact thermodynamic analysis of systems with few and many particles in terms of a set of ther- 
f~ l ' modynamic measures is discussed. We found that an exact thermodynamic description is feasible 

O ' based in the adoption of an adequate set of measures and the search of the free energy dependence 

Q^ , on the adopted measure set. A relation similar to the Laplace equation for the fluid-vapor interface 

S ■ is obtained which express the equilibrium between magnitudes that in extended systems are inten- 

sive variables. This exact description is applied to study the thermodynamic behavior of the two 
"j5 , Hard Spheres in a Hard Wall Pore for the analyzed different geometries. We obtain analytically the 

external work, the pressure on the wall, the pressure in the homogeneous zone, the wall-fluid surface 
tension, the line tension and other similar properties. 
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I. INTRODUCTION 

C \ The exact analytical evaluation of the partition function and thermodynamic properties in systems of confined 

particles is a new trend in statistical mechanics. Due to the inherent difficulties in searching the exact solution 
of three dimensional systems, the interest is focused in few confined particles, and is restricted to Hard Spherical 
particles. Systems composed by many Hard Spheres (HS) have attracted the interest of many people because of they 
constitute a prototypical three dimensional simple fluid |25| . Even, though its apparent simplicity only a few exact 
J^ , analytical results are known. In the limit of large homogeneous systems, only the first four virial coefficients in the 
_. ^ ■ pressure virial series for the monodisperse system are known (see ^28] and references therein). Similarly, the fourth 
r — \ order coefficient for the polydisperse systems were also obtained Q- It is interesting to note that the exact equation 
f^ of state (EOS) for the HS is unknown although an approximate, simple, analytical, and accurate EOS was found by 
Carnahan and Starling [ll|. The earlier published works on HS were specially devoted to the analysis of uniform 
fluid properties, as it was the classical Molecular Dynamic experiment on fluid particles of Alder and Wainwright (J|- 
^) ' Gradually, the focus of succeeding works turns to inhomogeneous systems. In the last decades a great effort were 
devoted to the understanding of HS inhomogeneous fluid systems, in part because such system are the starting point 
of several density functional theories |26l . |35| . These general theories deal with a large class of simple and complex fluid 
systems with successfully results in the study of the substrate-fluid behavior including wetting, capillary condensation, 
J^ ' and adsorption phenomena. Recent advances in the analysis of fluid adsorption in porous matrix were supported by 
H ' developments in this held [2y->Bl]. In last years much attention was focused to small systems of HS conflned in vessels. 
The study of simple fluids constrained to small cavities of various shapes has enlighte ning fundamental questions 
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of statistical mechanics and thermodynamics (for example about phase transitions |2ll . |31|), but only recently the 
relevance of few body systems was recognized. 

Few bodies conflned systems is a topic of statistical mechanics which belong at the opposite of the thermodynamic 
limit. The study of such systems is becoming technologically interesting because the manipulation of matter in the 
microscopic and nanoscopic scales shows that they can be built. Besides, the design of new nano-devices could take 
advantage of its properties. From that point of view, the use of simple hard-core potentials enables a schematic 
description of the interactions between particles and with the container. As we will see below, this simplifled picture 
makes analytically tractable the two- body sy stem. Interestingly, colloidal particles with HS-like interaction has been 
produced and studied experimentally |22l. l33l. |34| . 

In few bodies systems different ensembles are not equivalents. The correct ensemble to describe the properties 
of a given system is such that better simulates its real properties. Thus, the canonical partition function of the 



confined few-HS systems attempts to describe the statistical mechanics properties of this system kept at constant 
temperature. Besides, exact canonical ensemble studies of few bodies confined systems provides the building blocks 
for an exact grand canonical study of them. The grand ensemble is important because the statistical mechanics 
theory of macroscopic liquids is largely developed in such framework. We recognize that the absence of exact results 
for inhomogeneous fluids in this framework is an obstacle which difficult the theoretical improvement of the theory of 
liquids. Thus, we expect that in the near future the connection between exactly solved few bodies systems and the 
theory of macroscopic fluids can provide new theoretical insight. 

From now on we will focus on the analytical exact solution of few HS system in a pore making emphasis on canonical 
ensemble results. Until present only the two HS (2-HS) system was tackled. Recently, the canonical ensemble 2-HS 
confined in a spherical cavity was solved |43| . and also, the system confined in a spherical cavity with a hard internal 
core was evaluated [43 ■ In both works the principal result is the analytic expression of the configuration integral 
(CI), but the one body density distribution and pressure tensor were analyzed too. Studies of the same system in the 
framework of the microcanonical ensemble has also been done ^^. The present work (PW) is devoted to the exact 
solution of the statistical mechanic properties of 2-HS into hard wall simple pores in the framework of the canonical 
ensemble. We expose new results for the cuboidal, the cylindrical and the ellipsoidal cavities. We should mention 
that the microcanonical ensemble CI of 2-HS in a cuboidal cavity found in [42] is formally identical to that analyzed 
in PW for the same vessel. However, we present a different approach to the integral evaluation and a simpler and 
more explicit expression of the CI. We have checked that both solutions are equivalent. 

In section [H] we show how a hard wall cavity that contains 2-HS can be treated as another particle. There, we 
show explicit expressions of the canonical configuration integral for 2-HS into three pores of simple shape. We study 
the confinement in a cuboidal, cylindrical, and spheroidal, cavities. The obtained exact CI are functions of a set 
of parameters X which characterize the different shapes of the cavity. In Sec. IIIII we analyze both, the one body 
distribution function and the pressure tensor, for some of the studied cavities. In this Section we also obtain an analytic 
expression for the intersecting volume between a cuboid and a sphere which appears to be a novel geometrical result. 
Sec. llVl is devoted to the search of some universal features in the CI of the 2-HS system constrained to simple geometric 
cavities including the cuboidal, cylindrical, spherical, ellipsoidal, and also the spherical cavity with a concentric hard 
core. A discussion of how to obtain a thermodynamic description of the system by transforming the CI from Z2(X) 
to a more interesting description Z2(M) where M is a set of thermodynamic measures is done in Sec. [V] There, we 
find the equations of state of the 2-HS system in the studied cavities and obtain some exact results for the many HS 
system in contact with curved walls. Final remarks are shown in Sec. I VII 

II. TWO BODIES IN A PORE 

The canonical partition function of two distinguishable particles in a pore is Q2 — hr^Z2 being A the thermal de 
Broglie wavelength and Z^ the CI, which may be expressed as a three nodes graph 
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Here e(vi) = ei — Exp{—/3U{ri)) with i = 1,2, e(ri2) — ei2 — Exp{—(3ip{r 12)), ri2 = ri — r2, U is the external 
potential acting on each particle, ip is the interparticle potential, and the integration must be performed over the 
infinite space. The accessible region of space for the ith-particle, fi, is the region where e^ > 0, and its boundary 
is dft. In PW we assume that il and dil are the same for all (two) particles. The labeled P node in Eq. ([T]) that 
represents the pore is linked to the particles by the e^ bonds drawn with continuous lines in Eq. ([T]). Particles are linked 
each to other by the Cij bond drawn with dashed line. Pores with hard walls have e^ = {1 if r^ S il, and otherwise} 
and then the e^ bonds fulfils the in-pore condition. For hard spherical particles dj = Q{rij — a), where 8 is the 
Heaviside function, Vij — \rij\ and a is the hard repulsion distance (it is also the diameter of one HS). Therefore, the 
eij bond fulfils the non-overlap between particles condition being null if particles overlap each other. Both conditions 
are mandatory for the non-null value of the integrand in Eq. ([1]) . It is clear that Z2 for a 2-HS system confined in a 
hard wall cavity is by its nature a geometrical magnitude. This means that Z2 depends on a and a set of parameters 
which characterize the shape and size of the cavity. Therefore, Z2 is a piece of the bridge that links geometry and 
thermodynamics. We will return to this point in Sec. |Vl Before the evaluation of the integral ([1]) we may perform 
some simple Mayer type transformations on it. Using the general identity e — 1 + f we may replace the ei2 bond 
and/or the e^ bonds. The introduced /12 function is non-null only if both particles are overlapping while fi is null if 
j-particle is in the pore. We will draw the functions e^ and /12 with continuous line while we will draw the functions 



fi and ei2 with dashed hne. Following this procedure we obtain 




(2) 

where each graph with an articulation node can be factorized and easily evaluated [18j taking into account some trivial 
identities 

® O = Zi , (3) 

Cy<D = -2b, = -^a\ (4) 

Here Zi is the CI of the one particle system, —62 is the usual second virial coefficient, and a plays the role of exclusion 
radius. Note that Zi depends both, on the shape of the empty cavity and the HS size cr. In the first row of Eq. ([2]) it 
was explicitly separated the independent particle term Zf from the second term that concentrates all corrections to 
this simple picture. This term is 2Zib2{pore) being b2{pore) the first cluster integral with the complete dependence 
on the size and shape of the pore [18]. Therefore the first row of Eq. is 

Z2 = Zf - 2Zib2{pore) . (5) 

From an opposite point of view we may regard the p-node as it was a particle. This allow us to recognize that the 
right hand side term in the first row of Eq. ^ is part of the third virial coefficient of a fluid mixture |43l. |45|. In the 
second row of Eq. © was also extracted the first non- ideal gas term, — 2Z162, which contains the usual second virial 
coefficient for homogeneous systems. Therefore, the third term contains the nontrivial core of the problem involving a 
complex dependence on the pore's shape parameters. It hides the inhomogeneous system dependencies and takes the 
control over the entire density regime, from low density (or large pore size) to the close packing condition. Moreover, 
this term produces ergodic-non-ergodic transitions and dimensional crossovers. To make a contribution to the last 
graph in the second row Eq. ^, one particle must be outside of il (or inside Cl, the complement of ft) while the other 
particle must be inside of 51, and also both particles must be near each other. This explains that for large pores the 
term scales with the surface area of the container which is a measure of the size of dfl. Even more interesting, this 
graph remains unmodified if we turn to the conjugate system of 2-HS confined in 17, i.e. the graph is symmetric with 
respect to the in-out inversion. More explicitly, we introduce a partition of the euclidean space E'^ = f2 U J7 being 
Voo = Zi{Q) + Zi{rt) the volume of the space. The Eq. ([S]) is valid for Z, = Z2{^) and Zi = Zi{il.) as was already 
stated, but also for Z2 = Z2{^) and Zi = Zi(J7). This is the in-out symmetry of the 2-IIS system conffired in a hard 
wall cavity |43| . 

Now we concentrate in the evaluation of Eq. ([I}. In principle, the integration is over the positions of both particles 
(with a fixed pore position), however, it can be rewritten as an integration over the coordinates of the pore and one 
particle (by fixating the second particle). Hence, we ffist fix both particles coordinates and integrate over the pore 
center position which allow that both particles be inside the cavity. The result of the integration is the volume W. To 
build the region with volume W we can follow a simple geometrical recipe. Choose one point of the cavity, e.g. the 
center, and draw two cavities centered at particle-1 and particle-2 positions. The cavity-2 must be the translation in 
r of the cavity-1, i.e., they must be equally oriented. The overlap between cavities -1 and -2 is the available region for 
the pore center. In Fig. [T]we show a schematic picture for a cuboidal pore. The result of the first integration is the 
overlap volume W, the grey region (color online) defined by the overlap of cavities -1 and -2. At a second stage, we 
should integrate over the position of particle-2 with coordinate r. The integration domain is the region outside the 
exclusion sphere (ES) with radius a and inside the zone of vanishing W which defines the external boundary (EB). 
The EB is determined by the region enclosed by all the positions of particle-2 when we support cavity-2 on cavity-1 
and translate it in all possible directions by keeping in touch the boundary of both cavities. Through this padding 
procedure the obtained EB is the region enclosed by the dashed line in Fig. [T] The CI of the system reads 
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Figure 1: (color online) Representation of the integration procedure utilized in the evaluation of Z-z. The drawn points indicate 
particles 1 and 2. The solid line draw cavity-1 and -2, while the overlap region with volume W is shaded. Dashed line plots 
the ES and EB. 

where e(r) = Q(r ~ a). By integrating only the pore center position we find an unnormalized two body density 
distribution, g{v) — Z2^W(y)e{r^. Interestingly, Z2 and Wiv) of the 2-HS confined system are related to the CI 
of other systems as it is the confined stick-particle (or dumbbell) [4^ which may be obtained by a sticky-bond 
transformation. The one body density distribution p(r) will be analyzed in Sec. IIIII A simple consequence of Eqs. 
([Il [S]) is that Z2 depends on the X parameters (introduced by the e^ bond and W volume) , which characterize the 
geometry of each pore. Now we are ready to solve Eq. ([5]) for some simple cavities. As we mentioned above PW is 
mainly devoted to study the 2-HS confined system of distinguishable particles. Even so, at the end of Sec. |TT]we make 
a brief comment about the 2-HS system of indistinguishable particles. 

A. CI of 2-HS in a cuboid pore 

The empty cavity is characterized by the length parameters L^, L' and L'^. We introduce the effective cavity length 
parameters Li — L\~ a, which characterize the available space for the center of one particle and the dimensionless 
lengths li = Lil a with i ~ x, y, z. Then we obtain for the cuboid shaped pore Zi = L^LyLz, and 
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where r = {x,y,z) (see Fig. [T]). The EB is a cuboid with doubled length sides and the W(r) dependence turns 
convenient to integrate Eq. ^ over < a; < L^, < y < Ly, < z < L^ and multiply by 8. Although, W{r) is 
positive defined for any r and it is a non analytic function. For practical purposes we will extend analytically W 
to enlarge the integration domain outside the EB box to x > 0, y > 0, and z > 0. Assuming L^ < Ly < L^ it 
is necessary to analyze the integral ^ for different pore size domains in the parameter space X — {L^, Ly, L^) or 
the similar X = {lx,ly,h)- Introducing the directions: xy ex. L^x + Lyij, xz ex L^x + L^z, yz ex Lyij -\- LzZ, and 
xyz ex LxX + Lyjj + L^z we may distinguish eight different regions of X. 

Region 1 

The large pore domain is defined by the condition that ES is completely enclosed into the EB, i.e., that li> \ with 
i — X, y, z. The integral ^ splits into the simpler ones. 
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per py/a^—x^ px/cr'^—x^—y^ 

I2 = 8 dx dy Wdz , 

Jo Jo Jo 

— "5""^ LxLyLz — —(T (L^Ly + LxLz -f LyL^) 

8 1 

—a^Lx + Ly + L,)^-a\ 



(9) 



which in terms of dimensionless length variables is 



-—(7 ixf^y^z "o"*^ Y'X^y ~r f'x^z \ '"yf'z) i 

^a''{l^.+ly + h)~^a'' . (10) 



Z2 = Ii - 12 = ^? - 2 Zi62(cu6) . (11) 



Then, we find 



Last expression is similar to Eq. ([5]). 

Region 2 

The ii'S exceeds only two faces of the EB domain. Here, the exclusion sphere showed in Fig. [ijshould extend beyond 
the EB at most in one direction normal to the faces of the box. As far as this direction was labeled as x, then we have 
^x < 1, ^a > 1, and Iz > 1. We define the auxiliary integral la^, its integration domain is the spherical cup outside 
the EB box in the x direction 



T-Zx — S dx dy W dz , 

J L^ Jo Jq 



o" |-^ (1 - Ixf (5 + Mx + ^x) - ^ h^z (1 - ^^f (Ix + 3) + 
\ V^~^x {8 + 911 -2 It)- Ix (I - arcsin(;,)) 



(ly+lz) 



(12) 



Z2=Il-l2+l3x, (13) 

In the same sense, we define I^i with i = x,y,z. The integration domain of X^i corresponds to the spherical cup 
outside the EB in the i direction, which completes the description of the set of functions {T^x, ^3yi ^3z}- 

Region 3 

We consider the situation when ES exceeds only four faces of the EB, where the exclusion sphere must extend 
beyond the cuboidal EB in x and y directions but not in directions z and xy. In this case Ix l£ ^, ly l£ ^, h ^ ^ and 
ll + ll> 1. The CI is 

Z2=Il-l2+l3x+l3y (14) 

Region 4 

The next domain to consider is when ES exceeds all six faces of EB but any more. It goes beyond the EB in {x, y, z} 
directions but not in {xy,xz,yz}. Then, we need ^i < 1 and If + I'j > 1, for i,j — x, y, z with i ^ j, therefore 

Z2 = Ii — X2 + I^x + ^3y + ^3z ■ (15) 

Region 5 

In this region £'5' exceeds at four faces and four edges of EB. The sphere fall off the EB in {x, y, xy} directions but 
not in z. Then we have obtain l^ + ly < 1 and Zz > 1. We define the auxiliary integral T^xy, its integration domain is 



the right angle spherical wedge outside the EB in both x and y directions. Note that the edge of the spherical wedge 
does not cross the sphere center. In addition, we define X2xy i its integration domain is the space outer to ES and 
inner to EB. 



2:3.2; = 8 / dx dy Wdz, (16) 

JL^ J Ly Jo 



dx dy Wdz, (17) 



both integrals are related by 
For T2xy we found 



1-2xy —T-1~ 1-3x — 1-Zy + "^ixy ■ (18) 



1, 



^2xy --'{-^llll{^-ll-ll)+Y^l 



^h ^i-ll- II [8 + 9 (/2 + ID - 2 (?4 + ;4) + Qilf^ 

^ ^ _ +lz[n{lx,ly)+my,lx)]} , (19) 

n{u,v) = yi"^ (8 + 9m^-2m'') -uarcsin(u) + -'y(3 + 6u^-u'')arcsin( ^^^ ) . (20) 

15 3 Vvl-w^/ 

The straight forward generalization oi l2xy and T^xy in Eqs. (1161 I17p defines the set of functions {T2xy, X2xz, X2yz, 
I^xyi Izxzi ^3yz}- The zouc of the phase space with non-null integrand in Eq. ([1]), i.e. the available phase space of 
the system (APS) breaks or fragments in two equal unlinked zones because the pair of particles can not interchange 
its positions anymore. In this sense we refer to an ergodicity breaking in the canonical ensemble, which introduce an 
overall factor ^ = 1/2 in the CI, therefore 

Z2=^*{Il~l2xy) ■ (21) 

Note that ^ was not explicitly written in Eq. ^ and then a ^ = 1 value was there assumed. 

Region 6 

In this region ES exceeds at six faces and only four edges of EB. Here, the sphere should exceeds the EB in 
{x, y, z, xy} directions but not in {xz, yz}. Then, the region in the parameter space is Z^ + /^ < 1, Iz < 1, ^^ + '^ > 1, 
and ly + l"^ > 1. For this region the particles can not interchange its positions, thus, the APS breaks into two equal 
and unlinked zones. The ergodicity breaking introduce the overall factor ^ = 1/2, 

Z2=^*ill-l2xy+l3z) ■ (22) 

Region 7 

When ES exceeds at six faces and only eight edges but any vertex of EB we have the seventh region. Here, the 
sphere exceeds the EB in {x,y, z,xy,xz} directions but not in yz. The parameter domain is l^; < 1, Z^ + ^^ < 1, 
Ix'^^'z ^ 1> a^iid ly + ll > 1. Again, APS breaks but now into four equal and unlinked zones each one characterizing a 
set of microestates, which is non-symmetric under some of the symmetries of the cuboid cavity. This is a spontaneous 
symmetry breaking phenomena. The ergodicity breaking produces a factor ^ = 1/4, and the CI reads 

Z2 — £, * {II — X2 + I^x + 1-3y + 1-iz — 1-dxy — ^'ixz) , 

= C * (2^1 - 2:22;^ + Isz - 1-zxz) ■ (23) 



Region 8 

The last region considered is when ES exceeds at six faces, twelve edges but any vertex of EB. Then, the sphere 
exceeds the EB box in {x,y, z,xy,xz,yz} direction but not in xyz. Then, l^ + l^ < 1 for i,j = x, y, z {i ^ j), and 
l'^ + ly + ll > 1. With this conditions the APS breaks into eight equal and unlinked zones which also involves a 
spontaneous symmetry breaking. The factor introduced by the ergodicity breaking is ^ = 1/8, while CI is 

Z2 = £,* (2^1 — I2 + 2^3x + ^3y + 2^3z — ^3xy — ^3xz — ^3yz) ■ (24) 

Finally, in the case that ES exceeds the EB also in xyz direction, the partition function becomes null because both 
particles do not fit into the cavity. 

B. CI of 2-HS in a cylindrical pore 

Let us define the usual length parameters, height and radius, that characterize an empty cylindrical cavity ij^, R' . 
The effective cavity length parameters are then Lh = L'j^ — a, R~R' — a/2 and the dimensionless ones are given by 
h = Lh/cr^ R = s~^ = 2R/a. For the cylindrical shaped pore we have Zi — ttLji R^. As it was above mentioned, we 
need to know the volume defined by the intersection of two equal and parallel cylinders, W{r, Lfi,R). It is related to the 
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intersection of two disks of equal radii R and separated by a distance r, Wdisk{r, R) = 2R^ arccos(r) — r (l — r 
where r = r/{2R) by 

W{r, Lh,R) = \Lh-z\ Wdrskir, R) . (25) 

Note that VF is a well defined function of r only for the range < r < 1. The EB is a cylinder of double lengths and the 
W dependence turns convenient to integrate over < z < Lh, < r < 2R and multiply by An. The analytic extension 
of W for values z > will be considered when it becomes necessary. We need to analyze the integral considering the 
parameters X = (i?, Lh) which define the allowed pore size domain. Defining rz on Rf + {Lh/2) z we distinguish four 
regions. 

Region 1 

The large pore domain is defined by the condition that the ES is completely enclosed into the EB, i.e., that h> 1 
and R > 1. The CI splits into, 

Ii=47r/ rdr I W dz = Zf , (26) 
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—a V(2i?)2-CT2(_6i?4 + 5i?V2 + ^4) 



+7ri?2(2i?'* - 2R^a^ + cr"*) arcsin(s) , (27) 

^1 ( — ^ , i;3;a) is the Gauss hypergeometric function which can also be written in terms of complete elliptic 
[l!,H|- The CI is then 

Z2 = Ii - I2 . (28) 



Region 2 

The ES exceeds only the bases of EB, domain. Here, the exclusion sphere should go beyond the EB only in the z 
direction and then h < 1, R > 1. We define the auxiliary integral la^, its integration domain is the spherical cup 



outside the upper base of the EB 



I3;, = 47r / rdr W dz , 

Jq Jlh 



,S-2 



= i?^7r/45 |-256/is2(-2 + Ts^ + is^) E farcsin U yj\ - hA , < 

+15TT (3 - 12(1 + h^)s^ + 64:hs^) - 256/15^(1 + 2s^ - Ss'') F (arccos (h) , s^) 

+2s ^/{l-h^){l-{l-h^)s^) [-45 + 2(75 + 41h^)s^ - 24(-5 - 4h^ + h^)s^] 

+30(-3 + 12(1 + /i2)s2 + 8(-3 - 6/1^ + h^)s'^) arccos fs ^/l - hA \ , (29) 

where F{a, b) and E{a, b) are the mcomplete elliptic integrals of the first and second kind respectively |1]. The CI is 

^2=11-12+13. • (30) 

Region 3 

In this case ES exceeds only the curved lateral face of EB. The exclusion sphere should go beyond the cylindrical 
EB only in the f direction being h > 1 and R < 1. With these conditions particles can not interchange its positions 
producing that APS breaks into two equal and unlinked zones. We introduce the auxiliary integral l2r is given by 

l2r = 47r / rdr \ W dz , 

Jo Jq 

3 45 

[(32i?4 - 28i?V2 - 3(7^) E{s-^) + (-16i?4 - 8R^<j^ + 3a^) K{s-^j)] , (31) 

where K{a) and E{a) are the complete elliptic integrals of the first and second kind, respectively. We may also 
formally define I^r = ^2 ^ ^2r- The ergodicity breaking produces a ^ = 1/2 factor, being Z2 

Z2 = C * (II - l2r) - C * (2:1 - 12 + Isr) • (32) 

Region 4 

This region appears when ES exceeds both the bases and the curved lateral face but not the edges of EB. Therefore, 
the exclusion sphere exceeds EB in the {r,z} directions, but not in fz. In consequence h < 1 and R < 1, but 
/i^ + r2 > 1. As happens in Region 3, here the APS breaks into two equal and unlinked zones due to the ergodicity 
breaking, being ^ = 1/2 and 

Z2^^*{Il-l2+ Isz + Isr) ■ (33) 

Finally, if ES exceeds also in rz direction, the partition function becomes null because both particles can not fit into 
the pore. 

C. CI of 2-HS in a spheroidal pore 

The last CI that we evaluate in PW corresponds to the ellipsoidal pore. We restrict the study to cavities where 
only two principal radii are independent, i.e. to the revolution ellipsoids also called spheroids. Therefore, two distinct 
shapes the prolate and the oblate ones will be analyzed. Let us consider an effective cavity with spheroidal shape. 
The effective length parameters are the principal radii R and R^ where Re is the different radius. Dimensionless 
parameters are R = s^^ = 2R/a, C = 2Rc/(t, and A = Re/ R. For A < 1 we deal with the oblate, while for A > 1 
we deal with the prolate, spheroids. The configuration integral for one particle is Zi = ^R^Rc- The volume of 



intersection of two equally oriented spheroids, W{r, z), is related with the volume of the intersection of two spheres. 
In terms of WgphereiQ, R) = ^ i?"^(l — |r + ir^) with g the spherical radial coordinate, and r — g/{2R) we obtain 



Wir, z) = XWsphereWr-^ + {z/\)\ R) . (34) 

Function W is well defined in the domain < r^ + {z/XY < 4i?^. The EB is a spheroid with double length radii. 
For this pore geometry none analytic extension is suitable and the W dependence turns convenient to integrate over 
< z < 2i?c, < r < 2R, and multiply by An. We need to analyze the integral considering the allowed values of 
parameters X = (i?. Re) that define the pore size domain. We may distinguish three regions. 

Region 1 

In the large pore domain ES is completely enclosed into EB, i.e. C > 1 and R > 1. As it was above described, the 
integral splits into, 

Xi=47r/ dz Wrdr = Zf, (35) 



^0 



I2 ^ An dz Wrdr 

Jo Jo 
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which concerns to both, oblated and prolated spheroidal cavities. For A < 1, (—1 + A^)^'^arcsec(A) transforms to 
(1 — A^)^/^arcsech(A). Although, for A — > 1 we obtain (— 1 + A^)^/^arcsec(A) -> 1 which is consistent with the spherical 
pore result [43| . In terms of X\ and I2 we find 

Z2 = Ii - I2 . (37) 

Next regions concern the situation where ES exceeds EB. Under such condition it becomes necessary to make a 
separate analysis of prolate and oblate, ellipsoids. 

Region 2 (oblate) 

Here we consider an oblate ellipsoid, A < 1. In this region ES exceeds on top and down directions of EB, i.e. in the 
direction of the principal axis z, but not in f. Therefore, we consider C < 1, R > 1. We find that in this region it is 
simpler to deal directly with Z2, we obtain 



Z2 ^ Att dz Wrdr 

Jo JVcr'^-z'i 



144 Vl - A2 



Rv/R^~~T(3 + 16R2-4R'') -3(6R2-l)arcsech(R2) , (38) 



with Zmax = Aw ^ 1 In^"^ ■ We may also, formally define I2Z = Ii — Z2 and I-^z = I2 ^ 1-2z- In the case that ES also 
exceeds EB in the r direction the partition function becomes null because both particles can not fit into the pore. 

Region 3 (prolate) 

Here we restrict to a prolate ellipsoid, A > 1 . In this region ES exceeds in the lateral direction the surface of EB. 
Then, ES goes beyond EB only in f direction, but not in z, i.e. C > 1, R < 1. Under these conditions APS breaks 
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into two equal and unlinked zones and the ergodicity breaking produce the factor ^ — 1/2. Again, in this region is 
preferable to deal directly with Z2 



Z2 = C*47r/ rdr W dz , (39) 

where Vmax = \/<^'^ ^ {X2R)^ /y/X^ — 1. This integral was solved by splitting it in several parts, after some work we 
obtain 



Z2C = II - —TT^R^Xa^ + -^;-R 



9 2 24 V2 3A2 



R Vl - R2 (3 + 16R2 - 4R4 



144V-1 + A2 

3 6;^2^2 

+ , (6R2 - 1) [arcsec(A) - arccosfR)] , (40) 

144V-1 + A2 ^ Ji y J V yj , \ J 

Formally, we can define l2r —Ii — C^^-^2 and Xj,r='I-i-'I-ir ■ In addition, we note that if the exclusion sphere exceeds 
EB all around, particles do not fit into the cavity and then the CI becomes null. The Eq. PO)) is the last result about 
analytic expressions for the CI of 2-HS system confined in the studied cuboid, cylindrical and spheroidal, cavities. 

In PW we deal with a pair of distinguishable particles. Even sow, we make a brief discussion about the CI of 
a system of two indistinguishable HS (2i-HS). The canonical partition function of 2i-HS confined in a cuboidal, 
cylindrical, spheroidal and other shaped cavities are easily obtained from the CI of a 2-distinguishable-HS with the 
introduction of minor modifications. The first obvious change comes in the partition function definition because we 
must introduce the correct Boltzmann factor then Q2,ind — \-^^^Z2.ind- Secondly, we must analize the difference 
between Z2 and Z2^ind- In principle, expression ([T]) gives the starting point to define both, Z2 and Z2^ind- However, 
the evaluation of Z2 for the studied cavities involves the factor ^ that modifies Eq. ([T]) in some regions. We recognize 
that Z2,ind — Z2 in regions where no extra factor appears. A detailed inspection of the origin of ^ also shows other 
different situations. In some regions the ergodicity breaking appears because particles can not interchange their 
positions, but this makes non sense for indistinguishable particles. Therefore, regions where ^ = 1/2 corresponds to 
£.ind = 1. In other regions the ergodicity breaking also involves the spontaneous symmetry breaking, in this regions 
we find ^ind — 2^. In summary 

Z2,ind = £,ind Z2(^ == 1) , with ^„jrf = 2^ if ^ 7^ 1 , 

and S,ind = 2 if ^ = 1. Remarkably, partition function relates simply by Q2 md — Q2 if ^ 7^ 1- and Q2 md — jQs if 
5 = 1. 

III. LOCAL PROPERTIES: DENSITY DISTRIBUTION AND PRESSURE 

In principle, the partition function of the system holds in its global Statistical Mechanic properties. Such properties 
are presumably obtainable from some derivatives of Q2- This makes interesting the study of the analytical properties 
of Z2 which is done in Sec. IIVI Now, we are also interested in the local properties of the 2-HS confined system. 
Therefore, we study two functions, the one body density distribution p(r) and the pressure tensor P(r). We begin 
with a general brief description of the properties of p{r). For any pore shape, p{r) is [18, p. 180] 

p(r) = 2Z2^ei{r) /" e2(r2) ei2(|r - ral) dra , 

= 2^2-1 ei(r)(Zi-^2(r)) , (41) 

J2(r) = - f e2(r2) /i2(|r - r2|) dr2 , (42) 

where J72(r) is the overlap volume between the cavity and the ES (with a radius) at position r. This ES is produced 
by one HS-particle located there. The complete integral is j p{r)dr — 2. For an arbitrary r, J72(r) is positive and 
continuous but non-analytical and may be piecewise defined. When the particle is placed sufficiently deep inside the 
cavity all the ES is inner to the boundary. Therefore, for r such that the shortest distance to the boundary is greater 
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than cr, 1/2 (r) reaches its maximum value 1/2 (r) — 262- This means that for big enough cavities of any shape a plateau 
of constant density 

Po = 2^2-1 (Zi - 262) , (43) 

develops at a distance to the boundary grater than a. When r becomes nearer to the boundary the function J2(y) 
decreases and p{v) increases. For r outside of the cavity p(v) — 0, but we can define its continuous extension y(r) by 
dropping out the ei(r) term in Eqs. (|¥T|l . Outside of the cavity, for distances to the boundary greater than g^ y(r) 
becomes constant because J^iy) = 0. The J72(r) for the cuboid and cylindrical pores may be expressed by combining 
the 262 constant and the geometrical functions {i72a(r • a), J72afc(r • a, r • b), J72afcc(r • a, r • b,r • c)}, where \^a^h^c\ 
represent characteristic directions normal to the cavity boundary with inward normal versors {a, b, c}. The function 
J2a is the inner overlap volume defined by the ES and one boundary surface that intersects it, J2ah is the inner overlap 
volume defined by the sphere and two intersecting boundary surfaces, and J2ahc involves three mutually intersecting 
boundary surfaces. The short cut r-a and similar are the (minimum) distance between the ES center and a face of the 
boundary with normal inward versor a. Although, r • a extends to negative values when r is outside of the cavity. We 
may mention that inner overlap volume clearly identify a unique volume and then this description is non-ambiguous. 
When position r is on a cavity surface with simple curvature and away from other surfaces (a distance greater than 
c) 1/2 (r) — J2a(y • a = 0) which reduces to simple expressions. In such conditions, we have 1/2(0) = 62 for the planar 
surface, J2^sph^) = 62(1 — |s) for a concave spherical surface and J72,sp/i(0) = 62(1 + fs) for the convex one [43, 45|. 
For r on the lateral curved surface of a cylinder the analytic expression involving elliptic integrals is known [2^] . Its 
power series is, J2,cyi{0) = ^2(1 ~ |s — -^2^^) + 0{s^) and J2,cyi{0) = ^2(1 + f s + g^s"^) + O(s^) for the concave and 
convex cases, respectively. The question becomes even worse for the spheroidal pore surface, where we found analytic 
expressions of J72,sphrf(0) only for points on the poles and on the equatorial line. 

A. Density distribution in the cuboidal pore 

For the cuboid cavity the boundary surfaces are orthogonally intersecting planes. Therefore, in cuboidal cavities, 
J2a is the inner overlap volume defined by the ES and a plane that intersects it, J^2ab is the volume defined by the 
sphere and a right angle dihedron that intersects it, and J2ahc is the volume defined by the sphere and a right angle 
vertex. We must include a brief digression about the volume of intersection of a unit sphere and a set of mutually 
intersecting planes. As we are primarily interested in the cuboid we restrict ourselves to sets of mutually orthogonal 
planes with at most three planes. We introduce the function /Ca(r ■ a) which measures the volume of the spherical 
segment or spherical cap, defined by the intersection of the unit sphere at position r and a half-space with inward 
normal a. The vector r goes from a point in the plane to the sphere center. For a = x we have r • a = x with x > 
if the center of the sphere is in the positive half-space. For —1 < x < 0, 

dx / dy' / dz' , (44) 

X Jo JO 

= |(l-^x)2(2-x), (45) 

but Eq. P5|) is also valid in the extended domain — 1 < x < 1. Naturally, 

47r/3 = /Ca(r-a)+/Ca(-r-a), (46) 

/C,(-r-a)=/Ca(-r-a), (47) 

where the label a — —a corresponds to the inward direction —a. The function /Ca6(r • a, r • b) is the volume between 
the sphere and a right angle wedge when the edge cross the sphere. The wedge is defined by the quadrant determined 
by the intersection of half-spaces with inward directions a and b. As far as the center of the sphere does not lie on 
the edge this spherical wedge is different to the usual one. For a = x and b = y we obtain 



IC,y{x,y) = 2 / dx' dy' / dz' , (48) 



V 



2xy v 1 — x^ — j/2 — 2 arctan( 



xy 



+x (3 — X ) a,rccos( ) + y (3 — y ) arccos( — 



^1 - x2 - J/2 

(49) 
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where Eq. (P5|) applies for —l<x,y<0,x'^+y'^<l but Eq. (^^ is valid in the extended domain —1 < x,y < 1, 
x^ + y^ < 1. The half-length of the portion of the wedges edge inside the sphere is 1 — x^ — y^. The function 
/Cab(r • a, r • b) has the following properties 

K-air ■ a.) = K-abir ■ a,r -h) + K-abir ■ a,-r -h) , 

/Cfc(r • b) = JCabir • a, r • b) + ICab{-r • a, r • b) , (50) 

47r/3 = ICab{r • a, r • b) + K.ab{-^ • a, r • b) + lCab{Y • a, -r • b) + K.ab{-^ • a, -r • b) , (51) 

/Cab(r -a, -r -b) == K^i{v ■ k,~r -h) , 
Kab{-Y-k,v-h) = /Cab(-r -a, r -b) , 
/CQb(-r-a, -r -b) = /Cgb(-r ■ a, -r • b) , (52) 

being the Eq. ((5T|) a consequence of Eqs. (|46l [50 |) . The last /C function in which we are interested is /Ca6c(r-a, r-b, r-c), 
the volume defined by the sphere and a right angle vertex inner to the sphere, 

/•\/l~2''^~^'^ ^\/l-2:'^-z^ /.-^l-K'^-j/'S 

K.a;yz{x,y,z) = / dz' / dy' / dz' , (53) 

TT 1 

= g+xyz- - [/Ca;(a:;) +/Cy(y) +A:^(z)] , 

As happened before, Eq. ([55)) apply for — 1 < a;, y, z < 0, x^ + j/'^ + z^ < 1 but Eq. ([5^ is valid in the extended 
domain —1 < x,y, z < 1, x'^ + y^ + z^ < 1. Interestingly, we were unable to perform the direct integration expressed in 
Eq. (|53|l. although it was evaluated making a geometrical decomposition into simple terms. We obtain the following 
properties for /Ca&c(r ■ a, r • b, r • c) 

ICab{r • a, r • b) = /Ca6c(r • a, r • b, r • c) + ICabc{r • a, r • b, -r • c) , 
/Cac(i" • a, r • c) = /Ca6c(r • a, r • b, r • c) + ICabc{r • a, -r • b, r • c) , 
ICbc{r • b, r • c) = K.abc{r • a, r • b, r • c) + /Ca&c(-i" • a, r • b, r • c) , (55) 

47r/3 = /CQf,c(r • a, r • b, r • c) + /CQbc(-i" • a, -r • b, -r • c) + ICabc{-r • a, -r • b, r • c) , 
+/Ca6c(-i" • a, r • b, -r • c) + ICabc{r • a, -r • b, -r • c) , 
+ICabc{-r • a, r • b, r • c) + ICabc{r • a, -r • b, r • c) + JCabc{r • a, r • b, -r • c) . (56) 

ICabc{r • a, r • b, -r • c) = /Ca6e(r • a, r • b, -r • c) , 

/CQbc(r ■ a, -r • Id, -r • c) = A:„bg(r • a, -r • b, -r • c) , 

/Ca6c(-i"-a, -r -b, -r -c) = /Cabg(-r • a, -r • b, -r • c) , (57) 

where other identities similar to Eq. (|57p may be obtained by symmetry considerations. To the best of our knowledge 
the basic geometrical functions /Cab and ICabc are new results never published before. 

Using the /C functions we can complete the picture of p{r) for the cuboid cavity being that functions J72 and /C are 
related by 

J2a{r-a) = a^JCa{r-a) , 
J2ab(r- a, r -b) == a^JCab{r ■ a,r -h) , 
J2abc{r ■ a,r ■h,r ■ c) = a^JCabc{r ■ a, r ■ h,r ■ c) . (58) 
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with f = r/(T. We take the three orthogonal planes at 2; = 0, y = 0, and z = with inward directions x, y, and z, 
respectively. Thus, {x, y, z} represent the perpendicular distances to this set of planes. We assume a cuboidal pore 
such that Li > 2a (Region 1) and 0<a;<y<z<l, therefore 



J2{r) 



= < 



262, 

262 - ICx{-x) - ICy{-y) , 

262 - ICx{-x) - ICy{-y) ~ JCz{~z) , 

)<^xv{x,y): 

JCxy{x,y) - JCz{~z) , 

IC-j-ix) - ICxy{x, ~y) - lCxz{x, -z) , 

lCy{y) -K,xy{-x,y) - ICyAy,-z), 



x,y,z>l, 

a;<l, y>l, z>l, 

X, y < 1, a;^ + j/^, z > 1 , 

X, y, z < 1, x^ + y"^, x^ + z^, y"^ + z^ > I 

x^ +y2 < 1, z > 1, 

z, x^ + y^ < 1, x^ + z^, y^ + z^ > 1 , 

x^ + y^, x^ + z^ < 1, y^ + z^ > 1 , 

x^ + y^, y2 + z^ < 1, x^ + z^ > 1 , 



(59) 



ICx{x) - ICxyix, -y) - ICxyz{x, y, -z) , x^ + y^, x^ + z^, 
^xyz{x,y,z), x2+y2 + z2<l 



y2 + z2 < 1 , x2 + y2 _|_ ^2 > 1 



Following a similar procedure we can obtain 1/2 (r) for Regions from 2 to 8. In Fig. [51 we show three contour plot slices 
of p(y) for a cube with L = ha. From left to right of Fig. [2] the first slice shows the behavior of p(r) at half height of 
the cavity, the second one refer to a near wall position while the third one describe the behavior of p(r) on contact 
with the planar wall. The nearest line to the top-right corner of the slices corresponds to pa^ = 0.0159 , 0.016 and 
0.0162, respectively. The step in density between lines is Apa^ — 0.5 10""'. In Fig. [5] all the relevant characteristics 
of the density profile p{r) are apparent. We can observe the plateau of constant density at a distance a from the 
boundary and the increasing value of p{r) going from the plateau to the cuboidal cavity boundaries. Figure [3] shows 
a plot of p(r) for a given path in the same cubic cavity (L = 5a). There, the path is composed by several straight line 
parts. It starts at the cavity center (c), goes to the face (f) center, next to the middle of the edge (e), and next to 
the vertex (v) . The rest of the path follows other highly symmetric directions of the cube. We can observe here that 
even when p(r) is a piecewise defined function, it is continuous and also derivable (peaks appear because the path 
change its direction abruptly) . The minimum value corresponds to the plateau of constant density. For cavities with 
smaller size the extent of the plateau of constant density is more reduced. The effect of the higher confinement may 
be seen at Figs. |4]and[5l where the density distribution for a cubic pore with L = 2a is presented. From the left of 
Fig. [4] the first slice of p(r) is at half height of the cavity. Other two slices are similar to Fig. [2l The nearest line to 
the top-right corner corresponds to pa^ — 0.18 , 0.20 and 0.26, respectively. The step in density between lines is now 
Apa^ = 0.02. As can be seen in Fig. [S]the plateau disappears, because only for r at c the ES is completely inside 
of the cubic cavity. It is also apparent from a comparison with Fig. [3] From Figs. [2l [3l |4] and [5] we can also smell 
out the general behavior of p(r) for the 2-HS in cavities with different geometries and the effect of reducing the size 
of the cavity. 



B. Density distribution in the cylindrical cavity 

For the cylindrical pore the set of relevant functions are 

{J2z{r.z), J2f{r.r), :722r(r-z, r.f)} = a^lJC^ir.z), /Cp(f.r), /C2r(f-z, f-f)} , 



(60) 
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Figure 2: (Color online) Contour plot of the density distribution for a cubic pore with L = 5a. As gray becomes darker the 
density becomes higher. 
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Figure 3: Density distribution for a closed path in tlie cubic pore with L — 5a. 



where the cyhnder axis is in z direction and f is the radial polar versor. The inward normal to the lateral face is f = — f 
and r.r is the shortest distance from the sphere center to the lateral surface of the cylinder with radius R. Here, 
the functions {/Cz(r.z), ICf{r.r), /Czr(r-z,r.f)} are defined by translating to a cylindrical cavity the description made 
for the cuboidal cavity. The function ^^^(r.z) was already analyzed in Eqs. (|44j|47p . On the basis of the analytical 
expression for the overlap volume between a sphere and an infinite cylinder obtained in Ref. [23J (see Eq. (3) therein) 
we may obtain JCfiv.r) in terms of elliptic integrals. Some properties of these functions are 



47r/3 = /Cj=(r-f)+/C^(r -f), 



(61) 



JCz{r-z) = ICzf{r ■z,r ■r) + JCzr(r-z,r -r) , 
/Cf(r-f) == /Czr(r • z,r • f) +/Czr(-r • z,r • f) 



(62) 



47r/3 = /Czf (r • z, r • r) + /C^r(— r • z, r • r) + /Czr(r • z, r • f ) + JCzri^'r • z, r • f ) 



(63) 



We do not find an analytical expression for /Czr(r.z, r.r), which implies that we are not able to describe p(r) near the 
circular edges of the cylinder when (r • z)^ + (r • f)^ < 1. However, the exact value of p(r) on the edge is 



/C„.(0,0) = -/C^(0) 



(64) 



For the spheroid cavity we only found analytic expressions of J72a(r.a) for points on the polar axis and points on the 
equatorial plane, but they are not presented here. Functions J72f (r • f) and ICf{r ■ r) for the spherical cavity were 
obtained in [43], and for dimensions other than 3 in [43, 45j. These expressions enable to obtain p{r) near a concave 
or convex spherical surface. In addition, p{r) at the spherical pore with a hard core can also be obtained analytically 
using the same 1/2? (r • f) and /Cp(r • f). 

C. Pressure 

The analytic evaluation of the pressure tensor P(r), a symmetric tensor of rank two, is much more difficult than the 
evaluation of p(r) in an inhomogeneous fluid. For that reason we will not make a systematic search for each geometry 
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Figure 4: (Color online) Contour plot of the density distribution for a cubic pore with L = 2a. As gray becomes darker the 
density becomes higher. 
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Figure 5: Density distribution for a closed path in the cubic pore with L — 2a. 



confinement as was done in Sees. IIII Al and IIIIBI Even, we only make the complete evaluation for some simple cases. 
The relevant task of a detailed and systematic study of P(r) for 2-HS system near simple curved walls is planned 
to be presented anywhere. We focus on the evaluation of the pressure tensor P of Irving and Kirkwood [20] . The 
components of P for the 2-particle system are Pab (r) = P~'^SabP (r) + Pat (i"); with 



P^Ar) 



rl^F^^ / dt5{Y-ri+tYi2) 



(65) 



where r^ is the coordinate of the z-particle, ri2 
integration we obtain the identity 



ri ~Y2, r1 



ri2 • a, and F^^ = F12 • b = -^-^ -^- By direct 



dri2 ri2 



I {y,Yi,Y2) = / dt5{Y -Yi+tYi2) = r^2'^ ^5(fi2 - u)6 (ri2 - u) 



(66) 



with u = ri — r = uu. For a fixed r we introduce a set of cartesian and spherical coordinates with the usual 
convention for the polar angles i.e. rf2 = cos(6'i2) sin (^12) ri2, r\2 — sin (6*12) sin ((/)i2) ri2, and rf2 = cos(0i2)ri2. 
We can re- write Eq. ()65p . and for example, the P^^ component 



PP^. (r) 



e(ri)e(r2)(5(ri2 -a) ri2cos^(0i2) /(r,ri,r2) (fri(fr2 



(67) 



Using Eq. (p5)l . changing the integration variables to d'^u(Pr 12, expressing all the distances in a units and both 
variables in spherical coordinates i.e. d'^ri2 — ?'^2 sin(0i2) d?'i2 d(/)i2 rf0i2 and d^u = u'^ sin {(f)) du d(f> dd , and finally 
integrating on d^r 12, we obtain 



/3Pf. (r) 



Zo (7 



eAY — u e r - 



(1 - u) • u] cos^ (0) sin (0) 9 (1 - m) du d(j) dO . 



(68) 



Note that the range of u is 1. For r at a distance from the wall greater than 1 the integral /3P^ (r) becomes independent 
of r, because for all the available values of u in the integration domain we have e(r — u) = 1 and e[r — (cr — rt) • u] = 1. 
Therefore, for such r in the region of constant density (see Eq. (pS)) and comments therein) we find 



l3Pg (r) = Z2" V^ / cos^ (0) sin (0) 9 (1 - u) du d(f> dO , 



= Z^^2h2 



(69) 



The other components of the tensor are P^ = Py = Pj^ and P^ = Py^ = P^ = 0. This is expected because 
the pressure tensor in a region of constant density must be isotropic. The scalar pressure and the tensor relates by 
PP — (3tr (P) /3, where tr is the trace. Therefore, the scalar pressure in the region of constant density is 



/3Po =po + Z^^2b2 = Z^'2 (Zi - 62) . 



(70) 



A similar procedure was applied in |45J to the study of the 2-HS system in D dimensions. There, using a different 
definition of P^^, the authors obtained the same result for Pq. Pressure tensor near a planar wall can also be evaluated 
starting from Eq. (|68p . We consider a wall with inward normal z and an inner particle at a distance r.z = z with 
< z < 1. Integrating on a domain defined by |r — u| < 1, |r — (1 — u) • u| < 1, and < u < 1 we find the normal 
component 



/3P^ (r) = Z^^2T:a^ / cos^ (</)) sin ((/>) du 



Z:;'b2Z (3 - z' 



(71) 



16 




Figure 6: (color online) Position dependence of pressure tensor near a planar wall. We have drawn magnitudes related with 
scalar pressure P in continuous line, Pn in dashed line and Pt in dot-dashed line. The behavior of Pn {z) — Pt {z) (the 
integrand of Eq. (175^ 1 is shown in dot-dot-dashed line. More details about the plotted functions in the text. 



Such result can be easily checked. On one side, for an inhomogeneous fluid with planar symmetry we obtain PP^ (r) — 
Z^ 2 {Z\ — 362) which is independent of the position as it would be expected. On the other side, the fact that the 
contact value at the wall surface must be /3P^ {z — ^) — p(0) which implies P^ (z = 0) = 0. By following an identical 
procedure we find for both equal tangential components that 



PP^ (r) = Z2"Vcr3 / sin^ (</)) dud(i) , 



Z^\h^z [3- 



6 In (z)] 



For symmetry reasons the non-diagonal components are null. The scalar pressure near a planar wall is 

;3P(r) = p(z) + Z-^2h2Z [1 - In (z)] = Z-^2 [Zi - 63 (l + 2(1 - z^)l2 + z In (z))] . 



(72) 



(73) 



Finally, the wall-fluid surface tension of the 2-HS fluid in contact with a hard planar wall and the position of the 
surface of tension are 



/37 = /3a / [Pat (z) - Pt (z)] dz = -Z^^la^ 



(74) 



= —7 (T 



\Pn (2) - Pt {z)\ zdz^ -0.35 



(75) 



In Fig. [6] we plot together the position dependence for the pressure tensor components and other related magnitudes 
near a planar wall. The dependence with position is highlighted by plotting dimensionless magnitudes independent of 
Z2. We plot (/3Pfc(z) - p(z = 0)) Z2/2& with P^, = P, P/y, Pt and {PiPn{z) - I3Pt{z)) Z2/2& with continuous, dashed, 
dot-dashed and dot-dot-dashed lines, respectively. We see that at contact with the wall all functions go to zero with 
finite slope. For P, Pn and Pt the null value at z = is a consequence of the contact theorem. On the opposite, 
functions attain their definitive homogeneous value at distance a from the wall. Similar to the planar case, the 
spherical symmetry produce only two independent components P^ and P^ . We have obtained analytical expressions 
for the Irving-Kirkwood pressure tensor P near a spherical surface. This was done for convex and concave, surfaces. 
Even, the evaluation is not straightforward and therefore the study of the pressure tensor for the 2-HS system near a 
spherical wall will be presented in a future work. Near a cylindrical wall the components of P involve more complex 
integrals that we do not attempt to solve. 

Additionally, it is interesting to note a simple relation between pressure and density in the region of constant density. 
Recognizing that Zi plays the role of the system volume we can define the mean density p = 2/Zi. Therefore, from 
Eqs. P51 [75)1 we obtain the local compressibility factor in the region of constant density 



(3Po 
PQ 



1 



1 b 

2 p-^-b 



(76) 
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Figure 7: A simple picture representation of the Z2 X-space domain for all the studied pores. From top to bottom of each 
graph the volume decrease. 

This is a local EOS because describes the properties in certain location of the entire 2-HS system. In Sec. |V]we will 
study thermodynamic or global EOS. Expression (|76l) is very similar to the EOS of a (bulk) van der Waals system 
without the term of attractive force between particles. They differ in the 1/2 factor present on Eq. (1751) . which is 
related to the small number of particles of the 2-HS system. The Eq. ([75]) is valid for all the studied cavities, and it 
was also obtained for the equivalent system of confined 2-HS in dimensions D ^ 3. As it was suggested in Ref. [45j, 
it seems that Eq. ([75)1 is a universal feature of a 2-HS system confined in a cavity with Hard Walls of any shape and 
for all dimensions D > 1. We note that for a small enough cavity that produce a vanishing size density plateau the 
value of po depends on the geometry of the cavity. For a spherical cavity we have po — while in other cavities po 
assumes positive values. 



IV. ANALYTIC STRUCTURE OF CI 



The usual classical statistical mechanics links some global thermodynamic properties of any system of particles 
with some derivatives of ln(Z2), this idea will be discussed in detail in Sec. |Vl Now, we simple recognize that the 
analytical behavior of Z2 is related to the physical properties of the 2-HS. Therefore, the goal of this section is the 
study of the analytic structure of Z2 as a function of pore size parameters X, with the emphasis in the non-analytic 
domain. We are interested in investigate common features between cavities with different geometries. By including 
results from [43|, |45| we compare the CI for two hard spheres constrained by five different simple geometries: cuboid, 
sphere, sphere with a hard core, cylinder and spheroid shaped pores. A picture representing the structure of the 
domains for those Z2 is shown in Fig. [T] There, each box labeled with R (R-boxes) represents a region of parameter 
X domain studied in Sec. [H] as a separate case. The analytic domain of Z2 is the union of the (open) domains 
represented by the R-boxes. Straight line paths show the boundaries between adjacent zones, i.e. the non analytic 
domain of CI, while the broaden lines highlight paths of maximum symmetry (L^ = Ly = L^ for cuboid and Lh = 2R 
for cylinder). The stars distinguish the non- analytic domains involving the ergodic-non-ergodic transition. Dashed 
lines plot the crossover to systems with reduced dimension OD, ID or 2D; the 2D effective systems are represented 
with dark rounded-corner-boxes. The 2D limit for the spheroidal cavity has a different nature and we do not draw 
the box for this 2D limit. From Fig. [7] we can sort the structure of the Z2 analytic domains for the studied cavity 
geometries in an increasing order of complexity: sphere, spheroid, sphere+core, cylinder, and cuboid. The sphere 
is the simplest geometry, the cuboid results the most complex while the spheroid, sphere+core and cylinder have a 
similar degree of complexity. Moreover, if we restrict from the cuboid cavity to a cube, or from the cylinder to the 
sym,m,etric cylinder, its structure becomes much more simpler. This shows that the increment of the symmetry result 
in a decrement on the number of parameters in X. In summary, cavities with high (poor) symmetry and few (many) 
number of parameters X produce a simple (complex) structure. In Fig. [7] we identify several interesting common 
features concerning different shaped pores: (a) the large pore domain Rl, (b) its boundaries, (c) the Ri — >■ Rj, the 
signature of the ergodicity breaking, (d) the Rj — >■ 2D limit that exist in cuboid, cylinder and sphere+core pores, (e) 
the structure Ri — > Rj — > ID limit, and (f) the structures Ri — > Rj — > OD limit, Ri — > Rj — > CD limit and particularly 



18 



the last sequence Rj -^ OD limit. We now analyze the relevant properties for each case. 

(a) The large pore domain Rl 

Firstly, we concentrate in large cavities. The different analyzed geometries show that the large pore domain is the 
easiest to integrate and frequently the CI has a simple functional dependence. From direct inspection (see Eqs. ((8]fTT|) 
and also, Refs. [43,|43|) we note that for cuboid, spherical and sphere+core cavities the CI is a polinomy, but a more 
complex analytic dependence appears for the cylindrical and spheroidal pores. A comparison with two dimensions 
shows that the CI of the system of two hard disks into a rectangular cavity is also a polinomy, although for a circular 
cavity it is not true. From all the available CI we observe that Zi 62 (pore) of Eq. ([5|) naturally decompose in a universal 
way showing a simple dependence on basic geometrical measures of the effective pore. In terms of the volume notion 
V = Zi we obtain. 



y 62 (pore) = Vb2 



02/1 + i2Le + C2.1 + C22-—; 



(77) 



The constant coefficients 62 (see Eq. ([4])) and a2 — a'^n/S are independent of thepore shape. a2 appears in the virial 
expansion of the fluid-substrate surface tension and adsorption (referred as W2 J3H5l,'37|) and particularly, for a HS fluid 
in contact with planar and spherical walls p, SH, |4^. Besides the volume, in Eq. ((77|) we introduce other geometrical 
characters of the effective cavity, the area of the boundary A and the total edges length Le. In table U we present a 
comparison of the set {£2', C2,i; 02,2} for all the studied pore shapes, where the dependence on edges length, surface 
curvature and edge curvature is traced. We note that Vb2{pore) in Eq. (|77p for cuboid, sphere and sphere+core 
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Table I: Coefficients of Vb2(pore), dependence on the cavity shape for the large pore region. 



shaped pores involves constant coefficients {^2; C2,i; 02,2}- The coefficient £2 that multiplies Le has a unique positive 
value having the opposite sign to the preceding area term. Naturally, the edges are the area boundaries. Then, we 
saw the Le term in Eq. (j77p as a correction to the previous one. We interpret i2{cub) and i2{cyl) coefficients as 
being originated in the right dihedral edge formed by the intersection of two smooth surfaces. The C2,i is in general a 
slowly varying function of adimensional parameters s — a/2R and A = Rc/R. It is constant for cuboid, spherical and 
sphere+core pores. The negative constant C2,i(cm6) has a sign opposite to the previous edges term. From that we 
consider it as an end-of-edge correction which corresponds to the eight right vertex of the cuboid. Then, seeking for 
each vertex contribution we may write C2,i{cub)/a^ = —8/96 and therefore each vertex produce —1/96. On the other 
hand C2,i{sph) and C2,i{sph + core) are positive, i.e. they have the sign opposite to C2,i(cwfe), and also, they are not 
corrections to an absent edge term. Therefore, their nature is different to that C2^i{cub). Coefficients C2^i{sph) and 
c-2,i{sph + core) are originated on the curvature of the surfaces and their sign is opposite to the previous area term 
which corrects. Therefore, the surface curvature should produce a negative value for C2,i for both, a cylindrical and 
spheroidal pores. We introduce now the usual surface curvature measures, normal curvature j and Gaussian curvature 
fc, which take the values {j = R~^, fc = 0} and {j = 2R~^, k = R~^} for a cylinder and a sphere, respectively. We find 
thatc2,i{sph)/a^ = AR-"^ tt/IU ^ .JJ{R) 6'^^^ and C2A{sph + core)/ a^ = {.JJ{R) + JJiR-h)) S'^^'> ^ 2c2,iisph)/a^ with 

®7r [45]. For cylindrical cavities we find that 
is the curved lateral surface area, JJ{R) ~ 
Acurvj'^ = 27rL/i/i?, and for large radius C2.i(cy/) ~ AcurvR~'^- An unified description of cyl, sph and sph+core pores 
at large R is C2^i{cyl, sph, sph + core) a^^ — Acurv (f j^ + k) S^^' , but more complex dependence exist at C2,i(sphd). In 
fact, for large curvature radius and quasi spherical ellipsoids A ^ 1 we find C2^i{sphd) ~ C2^i{sphd) (1 + 4/5 (1 — A)^). 
Similarly, C2,2{cyl) relates with the curvature of the edges. We may resume some characteristics of {-F(s), G{s),H{X)}, 
F{s) and G(s) are positive and monotonically increasing functions in the domain [0, 1] with asymptotic minimum 
F(0) = G{0) = 1. H{X) is positive in its domain (0, co) and has a minimum at H{1) — 1. Its asymptotic behavior is 
H(X — >• cx)) — > A3Tr/16 and _ff(A — )• 0) — )■ A~^/4. In Fig. [8] we plot F{s) and G(s) adimensional functions. 

We have found a general structure of Vb2{pore) that is explainable by a hierarchy of correction terms. Term Vb2 
is the homogeneous component, it is linear in the volume and positive. The correction to Vb2 is the area term, the 
first signature of inhomogeneity. The area term is negative and then opposite in sign to the homogeneous term that 
corrects. Two types of essentially different corrections to the area term were found they comes from the edges and 
the curved area. The edge term which corrects the area term is negative and proportional to Le. For right dihedral 



the extensive quadratic curvature JJ{R) ^ Aj^ ~ 2'^'k and 6^^> ~ 3 ^2" 
C2,i{cyl)/a^ = A^urvR~^F{x)T:/?,M = JJ{R)6'-^^F{x)3/2, where A^,., 
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edges we found the value —1/15 for the constant of proportionahty, it appears for cuboidal and cyhndrical pores. The 
curved area term is a correction to the area term too, and sometimes, it is independent of pore size parameters being 
a constant. It is negative and approximately proportional to an extensive-like quadratic curvature AcurvP ■ This term 
appears at cylindrical, spherical, spherical+core and ellipsoidal pores. Noticeably, it does not exist any extensive-like 
linear curvature term. Two terms which correct the edge term were also found. They concerns an edge boundary 
term and an edge curvature one. Both of them basically reproduce the behavior of the corrections to the area term. 
These conclusions make interesting the evaluation of several coefficients in other geometric confinement, which may 
include, I2 for the edge of an arbitrary dihedral angle, C2.1 for a general vertex produced by three non-orthogonal 
surfaces, for the cone vertex, and the curvature correction of the general edge. 

(b) The boundary of the large pore domain Rl — > Ri 



In the rest of Sec. II VI our main purpose is to study the non-analytic behavior of CI when we go in the parameter space 
from an analytic domain to a contiguous one. With this in mind, we consider closed regions in the (Real) parameter 
space consisting in a region of the analytic domain with its boundary. We introduce the difference between the series 
representation of CIs, shortly 2'2(Ri) — Z2(Rj), corresponding to contiguous regions and evaluated in the neighboring 
of the common boundary. This may be not a well behaved magnitude. Even that, when at least one of the CI can 
be analytically extended in the contiguous region the difference Z2(Ri) — Z2(Rj) is easily analyzed. More complex is 
the case where neither Z2(Ri) nor Z2(Rj) can be analytically extended in the domain of the other. In such a case we 
made a careful comparison between the coefficients in each series. 

When we walk in the X-space from Rl to its outside the pore becomes unable to fit both particles for some fixed 
direction fi2. For example, going from Rl to R3 in the cylindrical pore becomes impossible that both particles locate 
in a plane orthogonal to the central axis. The effect on the volume of the available position phase space is not smooth 
enough producing the non analytic behavior of CI. We find that the behavior of CI in several paths of the type 
Rl — ;► Ri are well described by 



Z2(R1) - Z2(Ri) = AZ2 « -Ae^ 



(78) 



that is, for many situations we verify that CI has a discontinuous third derivative when the large pore domain is 
crossed in the parameter space. Here e = 1 — Li/ a^ is an adimensional vanishing parameter, e > and A/6 is 
the discontinuous step in the third derivative in the path Rl— >Ri. When ES exceeds the planar regions of EB, i.e. 
Rl— >R2, R1^R3 and Rl^-R4 for cuboidal pore; and Rl^-R2 for cylindrical pore, we obtain 



A(cm6) = {2TT(j^/'i)LyL^ , for Rl ^ R2 , 

= (27rcrV3) {LyL,^ + L^L^) , for Rl ^ R3 , 

= (2^aV3) [LyL, + L^L, + L,Ly) , for Rl 

A(cyZ) = (27rV^i?V3) , for Rl ^ R2 , 



R4, 



(79) 



where each equation should be evaluated at L^ — > cr consistent with the analyzed path. Here, the non-analyticity of 
Z2 is a consequence of the limiting behavior of the functions {I^^icuh), X^y^cub), I^zicub)} and I^^^cyl). Close to 



O 



„ 1.1 




Figure 8: Shape dependent coefficients in 62 (pore). Functions F{s) and G(s) of Table |l] are shown in continuous and dashed 
lines respectively. 



20 



the boundary they behave 



l3^{cub) w — naLyL, {a - L^f = -62 V^ (1 - L^Jaf , for Rl ^ R2 , 



Iszicyl) « -g'^^-R (^ - Lh) = -hV (1 - Lh/aT , for Rl ^ R2 , (80) 

where L^ ^ cr and i/i — > a for cuboidal and cyhndrical pores, respectively. The Eqs. (1791 I80p may be accomphshed 
with 

A = ^A,, (81) 

I3^ - -b2^ A+ (l - L,/af , (82) 

where A^ is the total area of such cavity boundaries which can not contain a sphere with a diameter. The same 
procedure is feasible for non planar boundaries, Rl— >-0D in the spherical pore, Rl— >-R2 in the sph+core pore, Rl— >-R2 
and Rl— >R3 in the spheroidal pore, and Rl— >R3 in cylindrical pore. Taking e = 1 — 2R/a we obtain 

4 
TTcr 
A(sp/i) = -—A , for Rl ^ OD , (83) 

which must be evaluated at i? = (t/2. The sph+core involves two non planar walls with different curvatures, the 
external spherical wall has radius R while the internal wall has radius i?i„. Both spherical walls are apart Lh = R—Rm- 
The gap in the third derivative with e = 1 — Lh/a is now 

4 
A(sp/i + core)\R, = ^ (A+ - 4^^) , for Rl ^ R2 , (84) 

where yl+ ~ Att{R^ + _Rf„) is the total area. 

We find three situations with a different behavior, they does not involve a finite discontinuity in the third derivative. 
The path Rl— )-R2 for the oblate-spheroidal pore has a discontinuous fourth derivative. For e = 1 — ^^ > we have 

AZ2 {sphrd) « ^^/(^''^ y2) s' ' fo^ Rl ^ R2 . (85) 

The path Rl — s- R3 involves an ergodicity breaking in prolate-spheroid and cylindrical, pores. Neglecting the factor 
^, for e = 1 — — > we obtain for the prolate-spheroid pore 

A^2 {sphrd) « ^'^' ^'''^''^' £7/2 f^^ Rl A R3 . (86) 

^ ' 105V-1 + A2 ^ ' 

We recognize that A{cyl) is somewhat ill defined cause their third lateral derivatives respect to 2R diverge logarith- 
mically to minus infinity. Even so, the difference between them becomes null. For £ = 1 — — >lwe obtain a 
non-analyticity expressible by the limiting behavior 

AZ^icyl) « _5^v^ ^7/2 ^ fo^ Rl A R3 . (87) 

Finally, the path Rl — > R4 in the cylindrical pore is analyzed by a superposition of results from Eqs. ([7^ [57)) . Its 
behavior is similar to that found in path Rl — > R2. 

(c) The path Ri — > Rj, a signature of the ergodicity breaking 

The rational power in Eqs. ([Ml [57)1 corresponds to path with ergodicity breaking, thus, we wish to study their 
characteristics. A third path with this behavior is R2 — > R4 for cylindrical pore. Again, neglecting the ^ = 1/2 factor 
we obtain the result described in Eq. (|57)) . based on the unanalicities oilsr- The cuboidal pore has also several paths 
of this type. They are the paths R3 A R5, R4 A R6, R4 A R7,R6 A R8, R4 A R8, and R7 A R8. All of them 
are characterized by the fact that a sphere with a radius fixed in the center of the cavity cross its edges. In fact. 
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this condition is equivalent to that described above for such a cavity (see Sec. Ill Al Region 5). Here the partition 
functions have an infinite discontinuous fifth derivatives as a consequence of the analytic behavior of the family of 

functions {Isxy, ^Sxz, ^3yz} 



^3xyicub) 



8 L, 
9!! LlLl 



j^~Ll-Ll)'^^ 



211/2 



aHa + a-^)^L. 



^l + ^l 



9/2 



for R3 A R5 , 



where a = L^/Ly and L+ is the total length of the crossed right edges, i.e. in Eq. ([55| L+ = AL^- Other paths are 
suitable analyzed by applying this result to the set {Isxy, ^Sxz, ^Syz}- The path R4 — > R6 is completely equivalent 
to R3 -^ R5. Somewhat different are the paths R4 — ?► R7, R6 — > RS, R7 — > R8, and R4 ^- R8, which involves an 
ergodicity breaking along with a spontaneous symmetry breaking. Even, their analytic behavior is basically described 
by Eq. §^. The path R7 A R8 is similar to R3 A R5 with the replacement y ^ z. Paths R4 A R7 and R6 A R8 
have two equal terms with the same value of a, the addition of both terms makes a unique contribution identical to 
Eq. ((88)) with L_)- the total length of the four crossed edges. Last path, R4 -> R8 involves three terms with a = 1, 
which resumes on one term with total edges length L+ = L — 12a/ V^. It is interesting to note that a similar situation 
is also possible for the cylindrical pore, where the circular edges are crossed by the sphere. It corresponds to the path 
R4 -^ OD which will be studied below. 

(d) The Rj -^ 2D limit 

The equivalent of the HS system in two dimensions is the Hard Disk system (HD). In the 2D-limit we may expect 
that 2-HS systems collapse to a 2-HD system. Then, Z2 should collapse to Z2,hd and then the CI of 2-HS in the 
cuboidal pore transforms to the CI of 2-HD in a box and so on. Expressions of ^230 for particles constrained in a 
rectangular or a circular pores, as well as, on the surface of a sphere are well known [27|, |4J, |4a|; this fact allow us 
check several results in PW. The expected limiting behavior of Z2 in terms of the vanishing length parameter e is 



Zo 



e' [Z2. 



HD 



^2,HDH?Ti 



-03+,(e) 



(89) 



where e — L^ and e — Lh for cuboidal and cylindrical cavities, respectively. Hence, we may study the unknown term 
£* Z2,}iDiim- For the planar surface 2D-limit we obtain q = 2, being for cuboidal shape 

1 



Z2.,WHrnicub) = - [nLyL, - 2a{Ly + L,) + a^] , for R2 ^ 2D , 



Lz 



a- -LI 



LyLz 'Arcsa\{Ly/a) - -L^ - aL^ 



for R5 ^ 2D , 



Ly v/ct^ -Ll+L, 



-'-^-^{-' 



Li + Li 



+LyLz (arcsin(i^/(7) + arcsin(Ly/cr) + 7r/2)] , for R7 -^ 2D , 



and for a cylindrical shape 
Z2,HV Hm{cyl) 



An 
135 



R^ 



1 



-1/2 



32R* - 157R2 - 3) + 45R2 arcsec (R) 



for R2 ^ 2D . 



(90) 



(91) 



In the case of a 2D-limit involving a curved surface confinement, we obtain for the spherical+core pore q — I and 

^23D ii,n{sph + core) = Att'^R [a^ - SR^] , for R2 ^ 2D , (92) 

where e = Lh- In the 2D limit of the oblate spheroidal pore R2 — >■ 2D we do not find the behavior depicted by Eq. 



(e) The Rj -^ ID limit 

The path going from Rl to the ID-limit has an ending structure Ri — > Rj — > ID. It means that, before to reach 
the limiting behavior a characteristic ergodic-non-ergodic transition appears. Once both particles are not able to 
interchange their positions the path Rj — ^ ID can happen and the final ID-limit may be attained. In that limit the 
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HS behaves like Hard Rods (HR) and Z-i collapses to ^2.hr- The limiting behavior for Zi written in terms of the 
vanishing length parameter e (e^ = L^Ly for a cuboid and e^ — ttR^ for a cylinder) is 

Z2 « e^ [Z2.KR + e" Z2MRHm] . (93) 

For the cuboidal pore ^2,hr = {Lz — c/2)^, q ^ 2, and 

Z2,}iRUm{cub) = —{Lz --a){a + a'^) , for R5 ^ ID , (94) 

OCT 

being a = Ly/Lx- For the cylindrical cavity we obtain Z2MR{cyl) = [Lu — cr/2)^, q = 2, and 

^2.HRi»m(c2/0 = — (Lft - a) , for R3 ^ ID . (95) 

TTfT 

In addition, we may compare with the ID-limit taken from the two dimensional 2-HD system confined into a rectangle, 
and from the 2-HD system confined between two concentric circles, from Refs. [23, |4J|. The ID- limit for the 2D 
rectangular confinement produces q = 2, while the circular pore with a hard core shows q = 1. We conclude that the 
power q = 2 is characteristic of straight line ID-limit while q = 1 corresponds to curved-closed-line ID-limit. The 
prolate spheroidal pore does not behave in accordance with Eq. (|93p . 

(f ) The Rj -^ OD limit 

The final state obtained in this limit consists of particles that cages in a final solid or densest configuration. This 
densest state of 2-HS characterizes by the complete spatial correlation of particles. Two different paths coming from 
Rl and ending at the OD-limit may be identified, they have the structures Ri — ?► Rj -^- OD and Ri — >■ Rj — >■ OD. The 
first case includes an ergodic-non-ergodic transition and sometimes also includes a symmetry breaking transition, it 
happens for the cuboid pore. We find that, in the OD-limit the phase space of positions (PSP) may collapse to three 
topologically different manifolds. For a cuboidal cavity the OD-limit shows a collapse of the PSP in a OD-manifold, 
i.e. a single point. Thus, the most compact state is a solid-like state. For the cylindrical cavity in the OD-limit 
the PSP collapse to a ID-manifold consisting in a simple closed line also called a circle. Here the densest state is a 
rigid body which is able to rotate with a fixed axis. For the spherical cavity the OD-limit shows that PSP collapse 
to a 2D-manifold given essentially by a spherical surface. Therefore the densest state behaves as a freely rotating 
rigid body. In the last two cases, even in the OD-limit, particles can interchange their positions. In general, the 
limiting behavior of Z2 in terms of some vanishing adimensional parameter e is Z2 (x s'' . For the cuboidal cavity with 
L = Lj. = Ly = a Lz and e = \/2 + a'^ L/a — 1 we obtain q — 6 and 

^^(-^) - ^^^-'^' + (2 + aW_i£!±0 ^a,. , f,, j,3 ^ 00 . (96) 

90a^ 105 

We note that q = 6 also in the case of a general cuboid. Analyzing the cylindrical geometry we find q = 9/2, 
L = Lh = a2R,e^ 2Ry/l + a^/a - 1, and 

Z2{cyl) = 7r^/2 (H-a^)^/^ 127575- V^e^/^ X [28350a-'^-^ , 

84105a-'' + 98100a-2 -f- 40200 -f- 6688a^ -I- 192a'*] , for R4 -^ OD . (97) 

For the spheroid, we can attain the OD limit in two different ways, by seeking the paths R2 -^ OD and R3 -^ OD. We 
obtain, q = 7/2, e = 2i?/cr - 1 > 0, and 

Z2isphd) « ^^^^^^ s^^ , for R2 ^ OD , (98) 

and also, q — 4, e — 2Rc/a — 1 > 0, and 

Z2 (sphd) « ^^^2^_i + ;^2) ^' ' fo'- R3 ^ OD . (99) 

The OD-limit in the spherical pore was previously studied in Ref. |43| . In that work, it was found q = 3 and 

£ — {2R/a— 1). Also, the OD-limit of a 2D system composed by 2-HD in a circular cavity has the same e but q = 5/2. 

We are now able to extract some minimal conclusions from this section. Based on the analysis made in (a) we note a 

very general decomposition of Vb2{pore) in terms of basic geometric magnitudes that characterize the effective cavity. 
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This decomposition could be applied in other confinement geometries. From (b) we find a common non-analytic 
behavior of Z^ when the ES exceeds planar regions of the EB boundary. It consists in a finite discontinuity at the 
third derivative with a step proportional to the surface area of the crossed planes. We also obtain a similar behavior 
for spherical surfaces and discontinuities at higher order derivatives in other curved surfaces. In general we observe 
that the paths between analytic domains involving ergodic-non-ergodic transitions Ri — ?> Rj are consistent with a CI, 
which scales with fractional powers of the vanishing magnitude. It is apparent in (b) where we find that a 7/2 power 
appears when ES exceeds a curved wall of the EB, and also, from (c) and (f) (see Eqs. ([551 Wi\ ) where we obtain a 
common non-analytic behavior of Zi when ES exceeds the right angle edges of the EB boundary given by a common 
power dependence of 9/2 in the vanishing length. 

A general picture of the dimensional crossovers agrees with the description given in [4q| . Given a iV-HS fiuid system 
in a region of the ZJ-dimensional space, the number of total spatial (i.e. translational) degrees of freedom is DF= N -D. 
When we consider a limiting process of dimensional crossover the dimension of the available space reduce to D ' with 
Q < D' < D. We define the number of lost degrees of freedom (LDF) as the power of the vanishing magnitude in 
the CI in the dimensional cross-over limit. We claim that LDF= N' ■ [D ~ D') where N ' is the number of particles 
constrained to the D ' dimensional region being usually N ' — N . One exception to this rule is the OD limit when the 
final densest state consists in a rotating N '-particle rigid-like system. In such a case we find LDF= N'D — n 3/2, with 
n indicating the number of independent degrees of rotational freedom for the caged N ' particles, being < n < I? 
[45| . In a unified description, for any dimensional cross-over we obtain 

LDF = iV' • (D -£>')-« 3/2, (100) 

where n = if D ' 7^ 0. Here, first term counts the lost of translational degrees of freedom while the second one 
compensates for the non- vanishing pure rotational degrees of freedom. For PW we must fix iV = TV' =2 with a 
starting value of I? = 3, and analize possible values D' =0,1,2. In the zero dimensional limit the 2-HS collapses to a 
dumbbell or stick. Thus, n = is a non-rotating stick, n = 1 corresponds to a rotating stick with fixed rotation axis, 
and ri = 2 is a freely rotating stick. Systems of two particles have a maximum value n = D — 1. Several sequences 
of dimensional crossovers described by Eq. (llOOp are accessible from the results exposed in PW. For example, in the 
cylindrical cavity the path R2 — ^ 2D involving LDF= 2 can be followed by a OD-limit with LDF= 5/2, obtained with 
£) = 2, £)' =0 andn = 1. 

V. THERMODYNAMIC PROPERTIES 

The aim of this section is to achieve the thermodynamic behavior of few bodies confined systems. Along this 
section we use the word thermodynamic in the sense of thermodynamic of fluids, where a fluid is a system of particles 
allowed to move in a given region of the continuous space. Our objective is to find the EOS that describe the global 
properties of a few body fluid system. In order to accomplish such a goal the discussion will be oriented towards the 
few and many HS system confined in a hard wall cavity with no restriction in the number of particles. In addition, 
we will keep in mind a system in a fiuid-like state. Besides these statements other systems could be included in the 
discussion without much effort, such as open systems and soft interactions. Again, we must emphasize that a few 
body system is far away from the thermodynamic limit N ^>- 00. Therefore, the thermodynamic description developed 
below does not concerns to such limit. In a few body system its different ensemble representations are not equivalent 
each other. Thus, we assume that the system under interest is well described by a certain Gibbsian ensemble and 
analyze the properties of this ensemble representation. From our point of view, we obtain the EOS of the system if we 
know the basic relations between the mean-ensemble values of the thermodynamic relevant magnitudes. A rigorous 
discussion about the equivalence between some mean-ensemble thermodynamic property e.g. U and the time average 
value Ur is out of the scope of PW. Even that, we can draw a general picture. We expect that for cavity's size in the 
ergodic regime and far from an ergodic-non-ergodic transition U = Ur for times r moderately short. For example, 
in a cylindrical pore it should apply in RI and R2, but far enough from R3 and R4 (see Fig. [5]). In case that the 
size of the cavity approaches an ergodic-non-ergodic transition the identity U = Ur only applies for increasing values 
of r. For cavities with sizes in the ergodicity breaking regime U and Ur may be different (e.g. R3 and R4 in the 
cylindrical cavity) . Next paragraphs are devoted to a general discussion about the thermodynamic description of few 
body systems, while at the end of this section we analyze the thermodynamic behavior of confined 2-HS systems in 
the canonical ensemble representation making a comparison between different shaped cavities. 

The pertinence of the thermodynamic theories to small systems was recognized by several authors, see e.g. the 
book of Hill p^] . From this book we can extract several arguments about the relevance of small systems to statistical 
mechanics and thermodynamics, and also, we find an interesting discussion about the particularities of the thermody- 
namics of small systems. Although, the central thesis of Hill is that the macroscopic thermodynamics must be adapted 
to extend its range of validity to include small systems. His thermodynamic approach begins with large (infinitely 
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extended) systems and drops to the small ones. Certainly, we adopt an opposite point of view. We state that the 
first law of thermodynamics concerns to few body systems, provided that, any assumption about the extensivity of 
the energy and entropy must be avoided. 

An implicit hypothesis of thermodynamics is that the equilibrium states of a large class of fluid systems may be 
specified with a unique small set of independent macroscopic quantities. A trivial example is the class of simple homo- 
geneous^ fluids usually studied by taking three independent macroscopic magnitudes (see e.g. Callen's thermodynamics 
book [fO] pp. 13 and 283). Therefore, we say that thermodynamics should have the Simplicity and Universality (SU) 
attributes. Usually, the studied systems involve a large number of particles, but does not exist a minimum cutoff in 
this quantity. To highlight this point, we note that in the statistical mechanics literature the grand canonical partition 
function is defined by a weighted sum of canonical partition functions over the available number of particles in the 
system (see e.g. [l^)- This sum starts from zero, following by one, two particles, and goes usually up to infinity. 
Therefore, systems with few bodies are included in the usual formulation of the statistical mechanics. We also note 
that usual relations that link statistical mechanics of partition functions and thermodynamic magnitudes do not make 
any assumption about the number of particles. This fact supports the idea that the same relations apply to systems 
with few bodies. Still, any assumption of extensivity in magnitudes like the energy, entropy, and free energies must be 
rejected in a few bodies system (see e.g. [10| pp. 360). We understand the thermodynamic pertinence of systems with 
many and few bodies as the Size Invariance (SI) of thermodynamics. Based on SU and SI, we argue that a consistent 
thermodynamic treatment of systems with large, many, and few number of particles should be possible using a basic 
small set of independent macroscopic quantities. Naturally, we will call to this the SUSI hypothesis. 

We want to bring attention to an unsolved problem in equilibrium Statistical Mechanics. At first sight it might 
be surprising that even when we may know the exact partition function of an inhomogeneous fluid system, their 
thermodynamic properties appear unrevealed. Our knowledge about the partition function comes from the exact 
evaluation of an integral (see paragraph above Eq. ([T])). As far as, the integrand and the limits of evaluation are 
functions of some set of independent parameters X, therefore by solving the integral we merely obtain Q(X). For a 
HS system in a hard wall cavity at constant temperature, the discussion is mainly focused on ZCX), where X can be 
of geometrical nature and usually involves proper lengths of the cavity, e.g. in a cuboidal pore X = {Ljc, Ly, L^}. 
Let us suppose that, for a given X space with dimension dimili.) the canonical partition function (5(X) for the N 
particles system is known within a reduced domain H. In such a domain we may obtain the Helmholtz free energy 

/3F(X) = -ln[0(X)], (101) 

which is related to other thermodynamic quantities by 

Ai = F(X) -i^_(X) , (102) 

F = U -TS , (103) 

dF =-SdT~dw . (104) 

In Eq. p02p the evaluation of the chemical potential /i assumes that the partition function for the system with N-1 
particles, (5_(X) (with Helmholtz free energy FJ) is also known in H. [/, 5, and T are the energy, entropy, and 
absolute temperature of the system, respectively. Lastly, dw is the differential of reversible work done by the system. 
Eq. (jl04p shows how F depends on both, T and X. The temperature dependence gives the entropy S 



dF 



(105) 

X 



while the X derivative at constant T is related to the work. Let us consider two different equilibrium states a and 5, 
characterized by parameters X^ and X(,, respectively. The variations Ai^, AS*, AC/ in going from state a to state h 
at fixed temperature are easily evaluated with the help of Eqs. p02l 11031 I105p . We may also evaluate the reversible 
work Wab hi going from a to 6 



Wab = - I VxF • dX = F{Xa) - F(Xfe) , (106) 

a 

where Vx is the gradient operator with respect to X parameters taken at constant T, and the line integral in Eq. 
(I106p does not depend on the path adopted between a and b. From here on. we implicitly make the same assumption 
for any derivative with respect to X. The Eq. (|106p enable us to define the differential of reversible work 

dw = dw^x = -%^ ■ dl = - Vxi^ • X d/ , (107) 
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being X some unit vector in the parameter space, dw^-^ the work to make a differential reversible change from Xa to 
Xf, = Xq + 'X.dl, and dj^ the directional derivative. Given any volume notion V, which may or may not be defined in 
the spirit of SUSI, we can define the overall pressure or pressure-for-work P^ j^ for an infinitesimal transformation of 
the cavity 

which makes sense only if VxV • X 7^ 0. For an infinitesimal transformation at constant volume we should ignore Eq. 
pOSp . Even, we may prefer to introduce some surface area notion A and therefore we can define an external surface 
tension or surface-tension-for-work by 

Eqs. (jlOSp or (|109p are indeed physical conventions, and therefore, we could describe the total work as it would be 
produced by either an effective pressure or a surface tension. From now on we assume that VxV • X 7^ 0. Then, the 
definition (|108|) is consistent with Eq. ()107|) . which now reads 

d^±-P^.,±dV^^ (110) 

where dVj^ = VxV • XdL The definition of P^ -^ requires the introduction of a volume notion V(X). Hence, pressure 
depends on both the adopted V and X. On the opposite, even when the choice of a different V modifies P^ -jj- it does 
not influence dw^-^. 

At this point we emphasize that, even when the above description is exact it is not completely satisfactory, ft 
says little about the thermodynamic properties of the fluid inside the cavity. It depends on X parameters, which do 
not have a universal thermodynamic meaning. The parameters needed to describe the shape of certain cavity are of 
different kind and quantity that those needed to describe other shapes. Even worst, for a given geometry they are non 
unique. We may extract some examples from the studied two particle systems. For a cavity with spherical symmetry 
we may utilize X = {R + cr/2} or X = {R}, but also, we may adopt X — {t^R^} all of them with dim{X.) = 1. In a 
cuboidal cavity we may adopt X — {L^, Ly, L^} or X — {Ix, ly, l^} with dim(X.) = 3, but also, if we are interested in 
a and b states with cubic symmetry we may choose X = {L} with dimpC) — 1. However, a somewhat more realistic 
cavity model may be adopted in which the substrate atoms, HS at fixed positions, are the building blocks of the rough 
confinement walls. In this case dim(X.) could be a much larger number. In addition, the X-representation prevents 
to compare results from dissimilar confinement conditions. Hence, the same fluid in a spherical or cuboidal cavity 
produces results which inhibit any comparison between them. 

We conclude that next step forward in the thermodynamic description of the system is out of the scope of the X- 
representation. Therefore, it is necessary to build the path between the X-representation of certain thermodynamic 
property, e.g. (5(X), and a universal description. Two basic questions have guided to us in the search of such a path; 
i) What properties of the confined systems should depend on the shape of the cavity? ii) What properties should 
depend on the particular choice of adopted parameters X? The rest of this section shows some answers, which arise 
from our inquiries. 

Being X an unsuitable set of parameters we must look for a better choice. At this point we wish to extract a 
paragraph from Callen's Themodynamics book, "It should perhaps be noted that the choice of the variables in terms 
of which a given problem is formulated, while a seemingly innocuous step, is often the most crucial step in the 
solution." ([l(l| p. 465). The interesting point is that Callen focus on the relevance of an adequate choice of variables. 
This question guide us to the concept of thermodynamic variable of state (VOS). We are interested in such VOS 
that characterize the spatial extension and other spatial features of an inhomogeneous fluid. A long time ago, in the 
origins of thermodynamics, volume was recognized as a good VOS for diluted gases as was stated in Boile's law in 
1662. A step forward was the introduction of surface area and curvature as VOS, it is documented in the study of 
vapor-fluid spherical interfaces made in 1805 and 1806 by Young and Laplace [2J,|43|. Although, in 1875 Gibbs (l2 | 
extended the use of curvature measures as VOS when he analyzed non-spherical fluid-vapor and fluid-fluid interfaces. 
Gibbs, also suggested the use of the length of the three fluid interface line as VOS. This idea was further developed 
in 1977 by Boruvka and Newmann [8|, which also introduced the curvature of such line as VOS. These VOS were 
extensively applied to the thermodynamic analysis in a variety of macroscopic inhomogeneous fluid systems including 
liquid- vapor and liquid- liquid interfaces, and adsorption of fluids on solids in accordance with SU |7l. Il4l Il6l Il7l. |39|. 
but they were never applied to the thermodynamic analysis of few body systems, in contradiction to SI. Besides, 
these thermodynamic magnitudes are based in geometrical concepts, but even when the geometrical concepts have 
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a precise definition, their counterpart thermodynamic magnitudes have usually not a precise meaning. For example, 
in the system of many hard spheres in contact with a (convex) spherical wall different choices for the locus of the so 
called Gibbs dividing surface is not innocuous. A comparison between Refs. Q and [3| shows that the locus of this 
surface may modify the volume and surface area of the inhomogeneous non-planar fluid system. Both modifications 
influence the macroscopic description of the entire system, changing the Laplace equation, the surface tension, etc. 
The most dramatic change is probably in the Tolman length. 

Therefore, we introduce a set M of thermodynamic measures, which should be suitable VOS in accordance with SUSI 
requirements. We seek for a set M with a precise definition which enables an exact description of few body exactly 
solved systems, and also, we expect that a good choice for M provides consistence with previous well stablished 
known results. The homogeneous fiuids are typically described by taking M = {V} with dim(M) = 1, while for 
inhomogeneous systems several authors currently add the surface area, being M — {V, A} and dim(M.) — 2. The 
classical analysis of the ideal gas produce an elementary EOS, PV — NkT. Accordingly, M must include a volume 
measure V with a pressure provided by P = —dyFCM.) compatible with the known system pressure, yielding the 
expected behavior for non interacting particles. The same thought applies for the surface area of the substrate and 
the wall-fluid surface tension 7. The discussion about the choice of M will be completed later in PW. Now, assuming 
that we have adopted a set M and also that M(X) is given, we must implement the thermodynamic description of the 
system using these measures. With this purpose we need to relate the X-representation and the M-representation. 
We state that Wab must be independent of the adopted representation X or M, then we claim 

Wab{M) = Wat{X) , (111) 

d^aM = dw,± , (112) 

where we assume that Ma and Mf, are well defined quantities and also, that for all X G H must exist M(X). Hence, 
Eqs. (fT06l [TOTll transform to 

b b 

Wab^- VmF ■dM = - mdM, (113) 



dw^^ = -d^F ■dl = -m-Mdl, (114) 

where m = Vm-F, and for a given direction X in the parameters space M — VxM • X. Comparing Eq. (|107p with 
Eq. pii)) we find 






dMi 



dX, 



X, , (115) 



where Xi is the i-component of X, rrij = -^^ , and M — M, means that all the measures but the I'-component 

are kept constants in the partial derivative. The Eq. (|115p is simply the chain rule for the F derivatives. When we 
adopt the volume notion of Eq. (|109l) as the volume measure Mi =V — V vje obtain P = — rni, and also from Eqs. 
(fTTClfTUgi) 



dim(M) 



VxMj • X 



P-P^,i, = ^P±= E ^^- S J ^ ^ (116) 



i=2 



VxV^-X 



is irrelevant 



which is a Laplace- like equation for a fluid-substrate interface [7J. The Eqs. (|109( I116P show that 

and therefore the restriction to unit modulus in Eq. (jllOl) is superfluous. An interesting point is that P^ ^ and P 
can be measured both experimentally and with molecular dynamic simulations. 

Now, to make a practical use of Eq. (|116p the unknowns nij, i.e. the EOS of the system, should be revealed. 
Therefore, we need i^(M) (see Eq. (|113p ). In general the set M may include dependent magnitudes and then 
d«77i(M) 7^ dimilC) showing that relation M 4-> X is not a one-to-one or biyective relation. Thus, the transformation 
F(X) — >■ F(M) is not a simple change of variables, which disable us to obtain F(M) = F(X(M)). We need a 
procedure to identify the hidden dependence of F{X.) in M. Accordingly, we must overcome two difficulties, find a 
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good set M(X) and obtain i^(M). Now, we can show that the selection of measures M and the identification of 
F(M) are not independent questions. To proceed, we analyze some results for the 2-HS confined system. 

We are mainly interested in fluid-like systems where particles can move freely and are able to interchange their 
positions. Then, we look for measures M that enable the thermodynamic description of systems in this regime. 
Certainly this M may or may not be suitable to describe other situations as solid-like or dense systems. The graph 
decomposition presnted in Sec. |ll](see also J4J|) and the analysis performed in Sec. llVIIVI show that some thermody- 
namic measures M appear naturally in F for cavity sizes in the Region 1. For higher confinement conditions, as in 
Region 2 and 3, some characteristics surface areas and lengths of the cavity also emerge as thermodynamic measure 
candidates. We focus on the results for Region 1 where any characteristic length of the cavity is greater than a. The 
list of measures candidates starts with the volume V = Zi — J e(r) dr, suggested by the graph decomposition in Eqs. 
([31 3]) . This volume appears usually in the study of inhomogeneous fluids [36-40] of different nature. Interestingly, for 
fluid systems in contact with hard walls, this V^(X) makes that P^ ^^ reduces to the contact pressure on the hard-wall. 
In fact, it reproduces exactly the hard-wall pressure contact theorem for planar, spherical and cylindrical hard walls, 
but also for much more complex geometrical shapes of the cavity [4J| ■ Other magnitudes are also suggested by the 
Eqs. (in [77l) . e.g. the surface area measure defined as ^ = J Weir) ■ fidr. We also consider Le, the measure of 
total edges length with right internal angle. More measures could be added, the number of right vertex, N^erti some 
measure of the surface curvature e.g. M = /(| j^ -I- k) dS, and a measure of the edge's curvature. 

Finally, even for Region 1, to ensure the exactness of Eq. (|116|) in principle we should include X in the set of 
measures. With all these measures we may conform a complete measure set Mc = {V, A,LeTNvert, M,X}, which is 
certainly not a small set of measures. We note that a hierarchy exist in Mc, the most important therm is V , the 
second in relevance is A. Both of them have been defined in detail, and its definition can be applied to a large class of 
systems. Next terms, Le, Ny^rt and M behave less important and their definition concern particular characteristics of 
confinement cavity. Finally, the last added terms to Mc are still less relevant. Their definition applies only to a given 
cavity geometry, and were included to make a complete description of F, so that Eqs. (|1111 I112p are guaranteed. The 
loss of relevance for incoming terms in Mc relates with the SU hypothesis. 

Now, we take into account all these questions to analyze the thermodynamic behavior of 2-IIS systems in Region 
1. The spheroidal cavity will be excluded from the thermodynamic analysis because we do not find a small set M 
that enable the unified study of this and other geometries. We adopt the small set of measures M = {V,A,Le,R} 
where the last parameter is the radius of curvature of the (curved) surface. Measure R is frequently used in the study 
of fiuid systems in contact with simple curved surfaces as such with cylindrical or spherical symmetries [7, 32]. We 
also select a rule to identify the dependence of F on the adopted set of measures M. It is based on re-writing Eq. 
(|77| in the form 

Vb2{pore) = Vb2 - Aa{R) + Lei{R) . (117) 

Again, the adopted M and the identification rule are non- unique. In the Appendix |X] a different M is analyzed. For 
the adopted M = {V, A, Le,i?}, we can define the F derivatives related with the volumetric- work, surface-area- work, 
edges-length-work and radius-of-curvature-work 



dF_ 
dV 



(118) 
T,M-y 



_ dF 

^ ~ 'dA 

- ^ 



(119) 



(120) 

TM-Le 



^ _dF 



(121) 
TM-R 



From Eq. (J116p we relate the difference of pressures AP-^ for an infinitesimal deformation in X direction with 7, r, 
etc. by 

^^x = 7 -^ + r -^ ^ + . . . , 122 
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Now it is apparent that Eq. (|122ll is a generalization of the Laplace equation obtained for a macroscopic fluid system 
in contact with a spherical wall ^14.. 15.]. An interesting fact is that the EOS given in Eqs. (|118m21|) may be strongly 
dependent on the details of the fluid system. On the other hand, the relation between AP-j^-, 7, r, etc. given by the 
Laplace-like equation (|122p only depends on the geometry of the cavity and the adopted M. For example, given a 
cuboidal pore it remains unperturbed if we confine 2-HS, an N-Lennard- Jones, or any other fluid. Before analyzing 
each confinement geometry we wish to state that, for all the studied cavities, the thermodynamic pressure from Eq. 

ma is 



f3P = Z^^2{V-b2) 



(123) 



This is our first global or thermodynamic EOS for the 2-HS system. The same expression was obtained in Eq. d75|) for 
the local pressure in the constant density region when we analyze cavities of any shape. We find that both pressures 
are equal, which shows the consistence of the present thermodynamic study. A similar result for spherical confinement 
was previously obtained [45]. Based on the universal behavior of Eqs. ([771lll7p and the consistence between the local 
pressure in the constant density region and thermodynamic pressure in all the studied cavity geometries, we confirm 
that the adopted volume measure is correct in the spirit of SUSL Therefore, taking the volume measure V — Zi we 
argue that the identity between both pressures should be true for a 2-HS system in any cavity shape. In the next 
paragraphs we perform the thermodynamic analysis for each pore shape. We fix a = 1 to keep notation simple. 



A. The cuboidal pore 



The cuboidal pore does not involve R, then M = { V, A, Le}. We obtain the thermodynamic pressure of Eq. (|123|) 
and also, 



/37 = -Z2-i2a2, 



(124) 



Pt^Z-^ 212 



(125) 



The three EOS relate the pressure, surface tension, and line-tension with the measures {V, A, Le} of the system. 
They apply to any cuboidal pore, in particular these equations are valid for the cubic confinement. Surface tension 7 
of Eq. (|124p is in coincidence with Eq. (fM)) . it is negative for large enough cavities. A Simple inspection shows that 
for large cuboids the EOS scales /3P ~ p -t- 62/3^/2, (S-f ~ —a2P^/2 and /3t ~ ^2P^/2 . For P^ -^ we may find in the 
literature two frequently used deformations. Adopting the length parameters X = {L^, Ly, Lz}, the first one is like 
a piston expansion transformation and reads X = (1, 0, 0). From Eq. (jlOOp 
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X 
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2{V- 



2a22{Ly + Lz)-2£2'i 



LyLz 



(126) 



with this choice of X magnitudes P, 7, r, and P^ -j^- are related each other by 



APv> 



Ly Liz 



LyLz 



(127) 



Equivalent results may be obtained with X = (0, 1,0) or X = (0,0, 1). The second option is an isotropic expansion 
with X = (1, 1, 1), which produces 
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(128) 



2Le 24 



(129) 



We now analize X = (1, 0, 0), X = (1, 1, 1) starting with a cubical cavity L = L^ = Ly ~ Lz- In this case Eqs. (|126p 
and p28p converge to a single expression. The same applies also to Eqs. (|1271 I129P which can be simplified because 
{A — Ax)/V — 4L~^, Lbx/V = 4L~^, and Le/A = 2L~^. Therefore, for all the studied X for a cubical cavity we 
obtain 



6P ^ = - 



Z-, 



2{V -b) + 



a2L- £2 ' 



(130) 
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APx-7^+r-^. (131) 

The same expressions ()1231 11241 11251 11301 I13ip are obtained if we start from the beginning the analysis of a cubical 
cavity with X — {L\ and X = (1), which shows the robustness of the procedure. 

B. The cylindrical pore 

From the same basic set of measures we recognize that the planar and curved surfaces, with areas A^ and A^,., 
respectively, are geometrically and therefore thermodynamically different. Then we split the total area in two, adopting 
M — {V, Ap,Ac, Le, R}. For P and jp we obtain expressions identical to Eqs. (|1231I124]) . Other EOS are 

/37, = -^2-^202 (1 " ciR-^) , (132) 

/3r = ^2-^24 (1 - ciiR"'^) , (133) 

PCr = -Z^^R-^ {cm Ac + civLe) , (134) 

where C/ = i^(s)/48, C// = G{s) 7r/14, cm = n/96 (F(s) + sF'(s)/2), and civ = 7r/210 (G(s) + sG"(s)/2). All of these 
coefficients are positive smooth functions with small values, e.g. cj < 0.03, which shows that 7p ~ 7,,. We can also 
extract the curvature dependence of 7 

--^^-^R-' + 0{R-'). (135) 

7p 48 

Taking X = {i/i/2, R} for this geometry, we see three simple choices for X. The piston expansion X = (1,0) provides 

^^x = > I • (136) 

For the lateral X = (0, 1) and isotropic X = (1,1) expansions we obtain 

iP„=,,|. + -,i + ,J^ + CH-i-, (137) 

AP^.,4+,.il±pi+.f + C«i. (138) 

For the cylinder with maximum area at fixed volume i/i/2 = R and Eq. (|138p reduces to 

APiM,, + 27.)A+(, + g)5|, (133) 

which may be still obtained analyzing the cylindrical pore from the beginning with Lh/2 = i?, X = {R} and X = (1). 

C. The spherical and spherical+core pores 

The spherical pore has Le — 0, which reduces the measures to M = {V, A, R}. The expression for P is identical 
to Eq. ()123p to which we add 

^7 = -Z2^2a2 (1 - 18"^i?"^) , (140) 

PCr^ -Z^^2a2A9-^R-^ . (141) 
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For X = {R} and X = (1) the difference of pressures is 



2 C, 



B. 



AP^-7^ + ^, (142) 

and the same result is obtained for any other choice of X and X. Note that using the Eqs. (|140[ I141[) we may 
transform ^ to obtain 

AP^=7(|+/3C.)+|I. (143) 

Eq. (|143|) is very similar to the Laplace equation for a macroscopic fluid in contact with a spherical wall p, [ij, [l^ ■ 
Now, we will analyze the 2-HS system in a spherical pore with a hard core. The shape of this pore involves two 
surfaces with different curvature. Therefore, as we did for the cylindrical cavity we consider two separate surface 
area measures. We adopt M — {V, A^, Ai, Re, —Ri}, where —Ri is the (negative) radius of curvature of the internal 
surface, and naturally, the labels e and i design the properties of the external and internal surfaces, respectively. The 
pressure P is given in Eq. p23p . while other EOS are 

/37e = -^2^^202 (1 - 18-^R7^) , (144) 

^7, = -Z^^2a2 (1 - 18-^R;^) , (145) 

/3Cfle = -Z^^2a2Ae9-^R~^ , (146) 

/3C_fl, = Z^^2a2A,9-^Rr^ . (147) 

In these Eqs. we recognize that the opposite sign in both radius of curvature do not affect surface tension expressions, 
but invert the sign in curvature term. Adopting the length parameters X = {Re, —Ri}, we can analyze three simple 
transformations X = (1,0), X = (0, 1), and X = (1, 1). For the first two cases we find 

(148) 



APx = 


2 , Crs 

-'^^Re+ Ae ' 


AP^ = 


2 ^ C-m 

^'R. + A. 



(149) 
ni J^i 

Here the effect of the negative curvature radius in the Laplace-like equations is apparent. The last transformation 
gives 

^ Re+R^ Lh Ae-A, ^ ' 

For all the studied simple geometrical confinement we obtained several relation between intensive-like magnitudes 
that resemble the Laplace equation. We must stress that previous to implement the thermodynamic study of the 
system we has needed to choose both, a set M and an identification rule (see Eq. (|117p ). Both choices affect its 
thermodynamic description. In the Appendix lAl a different choice for M is taken, which produce other set of EOS. 

In Figure [5] we plot pressure for work, pressure, surface tension and other EOS for the cube, the Lh = 2i? cylinder 
and spherical cavities (subfigures a, b, and c, respectively) as functions of the rough density. In continuous line we 
plot P^ -jj., in dashed line P, in dot-dashed line 7, dot-dot-dashed line shows r, and dotted line Ch- The vertical 
line shows the end of the plateau of constant density p^ and constant pressure Pq . The most remarkable feature is 
that Figs. [5^, [5)3, and [5}: are very similar in the density range (0, 0.5). A very small difference in P^ y^ at pa^ — 0.5 
is due to the expected geometric dependence of the ratio V/A. In Fig. [SJd both wall-fluid surface tensions, 7p and 
7c, are plotted but are indistinguishable. A clear difference between the three figures appears in the curvature term, 
which is not present for the cube cavity at Fig. |9^. In Fig. |91d Cr mix two curvature contributions, one due to the 
curved surface and other due to the curved edges. Even, in Fig. [S}: Cr is a purely surface curvature effect. A second 
difference is the vertical line that shows the end of the plateau, which is a purely local property of the 2-HS system. 
Note that P(X) (and also P(M)) is an analytic function at this point because their analytic domain extends to the 
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Figure 9: (color online) Magnitudes related to the thermodynamic work: pressures, wall-fluid surface tensions, line tensions 
and curvature-work. From top to bottom a) is for a cube cavity, b) is for a symmetric cylindrical cavity and c) is for a spherical 
cavity. The power n is chosen in each case to give a dimensionless magnitude. 

end of Region 1. The non-analytic point for Figs. [3^, [3]d andlHt corresponds to maximum densities pa^ = 2, 2.55 and 
3.82, respectively. Beyond the vertical plateau-end-line the identification of the thermodynamic pressure P with the 
plateau's pressure Pq breaks down because the central plateau of constant density disappears. If we wish to retain 
the identity beyond this point we can regard about the analytic continuation of Eqs. P5)) and ([70)) . This approach 
may conduce to a non-monotonic behavior of P related in some cases (for the spherical cavity) with a negative poi 
although, the total work P^ -^ is not influenced by this question. 

In consonance with Eq. (jl35p we visualize the possibility of analyze a cavity that mix planar and curved spherical 
surfaces. Such truncated-spherical cavity should have a Z2 involving a complex dependence on some set of parameters 
X. By virtue of Eq. ([77|) we infer that each surface makes its own contribution to Z2 allowing to obtain both the 
wall- fluid surface tension related to the spherical surface 7c and that corresponding to the planar one 7^. Then 7c 
should be essentially given by Eq. p40l) and 7,3 by Eq. (|124p with a common unknown function ^2, therefore 



^-1 
7p 



1^^" 



0{R- 



(151) 



Other interesting cavity is the half-cylinder, it mix curved and linear right anlgled edges. Even that we ignore Z2 
taking Eqs. (|125p and p33p we can obtain the curvature dependence of t 



— - 1 = R- 

T„ 14 



0{R- 



(152) 



The idea of build mixed shape cavities allow us to explore several confinement conditions involving complex geometrical 
shapes. As was already stated at [4^ and discussed in Sec. |lT]some results of PW are easily mapped from the 2-HS 
system confined in a bounded region to the 2-HS system confined to the conjugated unbounded region. This is a 
consequence of the inside-outside symmetry. Particularly, all the expressions which are independent of Zi and Z2 
are symmetric with respect to an inside-outside transformation. Therefore, Eqs. (|1351 11511) and (|152p may also be 
applied to the conjugated system where both particles are outside of fi. In general, the thermodynamic description 
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of the conjugated system is obtained by mapping Z\ and Z2 (see Sec. |TTl, and inverting the overall sign of P, P^ j^ 
and APjj.. 

D. Extrapolation to systems with many HS 

In [45| was recognized that some properties of the 2-HS systems can be mapped exactly to the many-HS systems 
in the low density regime. We simply follow the arguments of that work. In large inhomogeneous systems the 
thermodynamic limit is frequently considered, and sometimes, becomes convenient to introduce a mathematical surface 
where the surface tension is supposed to act. This is the so called Gibbs dividing surface. In our previous analysis we 
have not introduced a Gibbs dividing surface. Even so, if we are forced to define it we must assume that our Gibbs 
dividing surface is placed in coincidence with the surface of diverging external potential, e.g. for the spherical cavity 
it is the surface of a sphere with radius R. The wall-fluid surface tension of a HS fluid in contact with a curved wall 
and its limiting zero curvature value at the same density relate by 

^^ - 1 = -cvR-^ + 0{R-^) + 0(p) , (153) 

7c (00) 

where the geometric dependent coefficient is cv{sph) — 1/18 and cv{cyl) — 1/48 + 0{R^^). Both results apply to 
the HS fluid inside the cavity, but also, for the fluid outside the cavity. This symmetry is clear for the spherical 
surface, from the study of a fluid inside the spherical cavity with a central core, see Eqs. (|144l I145|l . Even, it is a 
consequence of the more general inside-outside symmetry. The central characteristics of Eq. (|153p is its zero order in 
density and second order in the radius of curvature R. We are now able to extract an interesting property of the HS 
fluid in contact with a curved hard wall. From Eq. ()153p the usual definition of the substrate-fiuid Tolman length 5 
(a magnitude independent of the radius of curvature) is 7c/7p — 1 = —2S/R + 02{R~^)- Therefore, we obtain 

5 = + O(p), (154) 

for both, spherical and cylindrical surfaces. It still applies for HS fiuid systems confined inside of the closed surface, 
and also, for fluids outside it. Our exact result for 5 is in contradiction with the constant value & = —a/A obtained 
in Eq. (35) of Ref. [9'| in the same limit. This difference would be consequence of the unusual volume definition 
adopted which does not reflect the volume available to the liquid's molecules (see Eq. (8) in ^] and the comment 
below Eq. (6) in Ref. Q)- In Ref. 01, Blokhuis et. al. analyze the behavior of a liquid system of particles interacting 
with a HS+attractive potential that mimics the London dispersion forces in contact with a curved hard wall. Using 
density functional calculations a limiting behavior of (5 ~ independent of the temperature is found (see Fig. 2 in 
that work) in good agreement with Eq. (I154p . In the same sense, our result for 6 agrees with Fig. 9 of Ref. |13J . 
Other magnitudes can also be evaluated. The line tension expressed to flrst non null order in density is 

/3r = p2|. (155) 

In consonance with Eq. (|152p . the flrst order correction on the line tension due to the curvature of the edge with a 
right dihedral angle is 

4^ - 1 = -^R-' + 0{R-') + 0{p) , (156) 

t(oo) 14 

which appears to be a novel result. The first non null curvature dependence for the density at contact is 

p{r = 0, i?) - p{r = 0, 00) ~ p^^cviR-^ ^ jf^^LLaR-^ + OiR-"") + 0{p) , (157) 

2 An 

with the packing fraction 77 = (t^/B) cr^p and cyi = 1,2 for cylinder and spherical cavities, respectively. For a convex 
wall we must invert the sign or simple change R — > —R. The Eq. (jl57p is in concordance with first density order of 
Eq. (36) in [9| which analize a fluid in contact with a convex hard cavity, but it is a new result for the HS fluid in a 
spherical cavity and also for the fluid in contact with a convex or concave cylindrical walls. 



33 
VI. FINAL REMARKS 

The analytical evaluation of the canonical partition function for the 2-HS confined system were performed for 
several cavities with simple geometry. The cavities considered were the cuboidal, cylindrical and ellipsoidal pores. 
The obtained expressions cover all density range from infinite dilution to the jammed densest configuration. The one 
body distribution function and pressure tensor were also analyzed. As a byproduct, we have obtained expressions 
for the volume of intersection between a sphere and a dihedron with right angle, between a sphere and a right-angle 
vertex, and thus the expression for the intersecting volume between a sphere and a box. To the best of our knowledge 
this expressions were not previously published. The three studied cavities were compared with the spherical and the 
spherical with a hard core, cavities, hence the study of simple pore's geometry is completed. The general behavior of 
all the available CI where analyzed by a graphical representation, which shows how the X-parameter space breaks in 
several open analytic domains. Attention was also paid to the CI solution for large cavities, to the characterization 
of the non analytic domain and the dimensional crossovers. 

Finally, we have focused on the thermodynamic properties of the 2-IIS confined system. Several questions about 
the free energy dependence on geometrical parameters X and its thermodynamic meaningfulness were discussed. We 
show the necessity of introduce a set of thermodynamic measures M based in extensive-like magnitudes. These neat 
defined measures constitute the basis of a consistent method developed to make the thermodynamic study. We find 
that pressures, surface tension and similar intensive-like magnitudes are then obtainable analytically. A common 
feature was the arising of an exact expression resembling the Laplace equation, which establishes the equilibrium 
between these quantities. Finally, several connections to the many-HS system in contact with curved hard walls were 
found. We have evaluated the first curvature corrections to the surface tension, Tolman length and line tension in the 
low density limit. 

The solved integrals to obtain the CI of 2-HS system also shows the complete dependence of b2{pore), the first 
non-trivial cluster integral, for the many HS system in the cavity but also outside it. For large enough cavity 62 (pore) 
is analytic. Even that, for smaller cavity's size h2{pore) is a non analytic function of the X and M-parameters. We are 
convinced that any cluster integral bj(pore) behave the same behavior. Cluster integrals are basic functions appearing 
in the virial expansion of the so called real gases EOS, thus, the study of the unanalicities of bj(pore) could be of 
interest. 

The performed study of free energy dependence of two-body simple systems on the geometry of the container does 
not close the prospection. Indeed, it shows that next steps should focus in the free energy contributions of dihedral 
edges (straight and curved ones), non right vertex and cone vertex. One of the conclusion of PW is that this future 
inspection should be numerical. 
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Appendix A: Thermodynamic properties with a different M 

In Sec. |V] we have discussed why the thermodynamic properties of the systems depends on the adopted set of 
measures. Here, we investigate on this dependence adopting a different set of measures to that used in Sees. IV Al 
IV Bl and IV CI Now, we utilize a measure set M currently used for the study of systems with spherical interfaces. We 
adopt M = (V, A), and add to this the parameters used for the description of the studied geometric confinement. As 
far as, we will restrict ourselves to the full symmetric cavities, cubic, cylindrical and spherical only one characteristic 
length parameter is needed. The unique parameter is X = {R}, where L = 2R for the cube and i/,, — 2R for the 
cylinder. Following the same ideas depicted in Sec. |Vl taking an area common factor in Eq. (|77p it transforms to 

Vb2{pore) = Vb2 - A {a2 - £2Le A'^ + C2,iA-^ - C2,2Le R-^A-^) ^Vb2-Aa{R) . (Al) 

Here, the expression between parenthesis only depends on R. We then analyze with M — {V, A, R) the isotropic 
expansion X — (1). For all the cavities we obtain Eq. (|123p for the pressure and 

(3j = -Z^'2a{R), (A2) 
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PCr^-Z-^2A^. (A3) 

We may highlight that a{R) is a known function and therefore 7 and Cr are analytically known in the three analyzed 
cavities. The Laplace- type equation is 

APx=7| + C«l (A4) 

Last expression without the jSCr term was obtained in some refined studies of spherical cavities in the bulk of fiuid 
systems, and also, in studies of spherical drops surrounded by its vapor. Notably, in our systems, which does not need 
to be spherical PCr ^ R^^ a higher order term in the characteristic length of the cavity. 
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